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Abstract 



We review an extension of Migdal's Theory of Finite Fermi Sys- 
tems which has been developed and apphed to cohective vibrations in 
closed shell nuclei in the past ten years. This microscopic approach is 
based on a consistent use of the Green function method. Here one con- 
siders in a consistent way more complex lplh(8)phonon configurations 
beyond the RPA correlations. Moreover, these configurations are not 
only included in the excited states but also explicitly in the ground 
states of nuclei. The method has been applied to the calculation of 
the strength distribution and transition densities of giant electric and 
magnetic resonances in stable and unstable magic nuclei. Using these 
microscopic transition densities, cross sections for inelastic electron 
and alpha scattering have been calculated and compared with the 
available experimental data. The method also allows one to extract 
in a consistent way the magnitude of the strength of the various mul- 
tipoles in the energy regions in which several multipoles overlap. We 
compare the microscopic transition densities, the strength distribu- 
tions and the various multipole strengths with their values extracted 
phenomenologically. 

PACS numbers: 21.60.-n, 24.10.Cn, 24.30.Cz. 
Keywords: Microscopic theory, Giant resonances. Unstable nuclei. 
Transition densities. 
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1 Introduction 



1.1 Aims of the review 

The atomic nucleus is a complicated quantum mechanical many-body system 
with a rich excitation spectrum. The experimental investigations of nuclei 
during the past 50 years have provided us with an overwhelming amount 
of excellent data which have to be understood in the framework of a quan- 
tum mechanical many-body theory. Moreover, we expect new data from the 
coming radioactive beam facilities for atomic nuclei far from the stability 
line, which will give us information about, for example, nuclei with a large 
neutron-to-proton ratio. Some of the features of the nuclei, such as their level 
fluctuations, can be described in terms of random matrix theory, which is now 
interpreted as the manifestation of chaotic motion. This stochastic nature of 
the levels is, however, only one aspect of the spectrum of the atomic nucleus 
which shows, on the other hand, many regular features, such as single-particle 
structure and collective modes. This regular behavior is a consequence of the 
Pauli principle and the special structure of the nucleon-nucleon interaction 
that leads to the concept of a mean field. As a consequence of the mean 
field, one is able to define quasiparticles in the sense of Landau and to apply 
Landau's theory of Fermi liquids P to the nuclear many body problem, as 
has been done by Migdal in his theory of finite Fermi systems (TFFS) j2]. 

In this article we focus on the question of how nuclei respond to a weak ex- 
ternal field. This corresponds in lowest order to linear response and Migdal's 
equations look formally very similar to it. As we will show in the next sec- 
tion, however, the range of validity is much larger than one would expect 
from the conventional derivation of the linear response equations. Migdal's 
equations are derived within the many body Green functions (GF) theory. 
Using this powerful formulation of the many-body problem one is able to 
obtain equations which in principle are "exact" . This is due to a renormal- 
ization procedure similar to the one used in quantum field theory. In the 
original version of the TFFS one considers explicitly only the propagation 
of particle-hole pairs in the nuclear medium. All the other configurations, 
such as the two particle-two holes ones, are renormalized into an effective 
two-body interaction and effective operators. The final equations have the 
same form as the corresponding equations of the conventional linear response 
theory, or random phase approximation (RPA). 

The important observation by Migdal was that the effective interaction 
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and effective operators so defined depend only weakly on the mass number 
and the energy. Therefore, as in Landau's original theory, these effective 
quantities can be parameterized and the corresponding parameters should 
be the same for all nuclei except for the lightest ones. Moreover, these pa- 
rameters have been chosen to be density-dependent, as the interaction inside 
and outside of nuclei may be quite different. These assumptions have been 
shown to be correct in numerous applications and they were and still are the 
most important features of TFFS. One can hope that the same parameters 
are also applicable to nuclei far from stability where we shall next apply the 
theory. 

It is well known that pairing correlations are important in non-magic 
nuclei. In this case on has to define quasiparticles analogously to the BCS 
theory, which include the pairing gap, and in the linear response equation one 
has to consider as well the change of the pairing gaps [21 E] in a consistent 
way. The equations of the extended theory look formally identical to the 
quasiparticle RPA (QRPA) g]. 

The standard TFFS is not a self- consistent theory in the sense that 
one starts with an effective interaction from which one obtains simultane- 
ously the mean-field quantities such as single-particle energies, single-particle 
wave functions and the particle-hole interaction for the linear response. In- 
stead, one has independent parameters for the single-particle model and the 
particle-hole interaction. Of the various extensions of Migdal's original the- 
ory that have been made over the years, the first were directed at taking 
self-consistency effects into account El El • 

It has been shown that the TFFS can also be applied to strongly de- 
formed nuclei [Sj. In that case the appropriate quasiparticles can be deduced 
from the deformed shell model of Nilsson. The extension of TFFS to a second 
order response theory has been developed and the formalism for the calcu- 
lation of electric and magnetic moments of excited states in even-even nuclei, 
transitions between excited states, as well as isomer shifts of rotational states 
in deformed nuclei, can be found in refs. fOl . 

In addition, the A-isobar has been included in the space of quasiparticles 
and, in the effective interaction, the effect of the one pion exchange, 
which turned out to be crucial for the existence of pion condensation and the 
properties of spin modes [71 [T21 (El . 

As mentioned earlier, both the linear response theory and (effectively) 
the TFFS are based on the RPA, which is restricted to a configuration space 
which includes one particle-one hole (Iplh) configurations only. If applied 
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to collective modes in nuclei it gives values of, for example, the centroid en- 
ergies and total strengths of giant resonances in good agreement with the 
data. The widths of the resonances and their fine structure, however, are 
not. These can only be obtained if one includes in the conventional RPA and 
its related approaches the coupling of the Iplh excitations to more complex 
configurations. Therefore, if one wants to apply TFFS to a realistic theo- 
retical interpretation of giant resonance experiments one has to extend the 
formalism to include explicitly in the theory more complex configurations 
than Iplh pairs. 

The main aim of this article is therefore to discuss such extensions of 
Migdal's theory and their application to giant resonances in nuclei. The 
essential feature — that which determines the necessity and the title of our 
approach — is of course the explicit taking of complex configurations into ac- 
count. After a short review of the standard TFFS we discuss in a qualitative 
way possible extensions of the theory. Here we will give physical arguments 
for why the coupling of the most collective low-energy phonons to the single- 
particle and single-hole propagators is the most important and which higher 
configurations are necessary for a quantitative description of the widths of 
giant resonances. 

In section 2 we derive the basic equations of the extended TFFS (ETFFS) 
for closed shell nuclei. As in the original TFFS, the GF method is used, and 
for convenience we restrict ourselves to doubly closed shell nuclei, where 
pairing correlations can be neglected. As we shall see, there exist several 
stages of sophistication of the extension, and we will discuss these in some 
detail. As a first approximation we include the complex configurations in 
the excited states only. This already gives rise to the fragmentation of the 
multipole strength, but the total strength remains unchanged. In the next 
step we include also the more complex configurations in the ground state of 
the nuclei. This gives a further fragmentation but, in addition, gives changes 
in the magnitudes of the electric and magnetic transition strengths, which 
turns out to be important for a quantitative comparison of the theoretical 
results with the data. The final formulas of the extended theory are given 
in r-space, which is especially appropriate for including effects of the single- 
particle continuum. 

In section 3 we apply the various versions of our theory to giant multipole 
resonances (GMR) in medium and heavy mass nuclei. The comparison with 
experiment demonstrates the power of our new theoretical framework and, 
simultaneously, the importance of the different steps of our approach. 
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In section 4 the calculations of transition densities within our approach are 
reviewed. With this information one is able to derive cross sections that can 
be compared directly with electron and alpha scattering experiments. One 
major difference between the microscopically calculated transition densities 
and the phenomenological ones is an energy dependence of the microscopic 
ones. In this connection we discuss also the problems that arise in the analysis 
of giant resonances in nuclei in which the resonances are very broad and where 
the various multipolarities overlap. 

Finally, we summarize our review and discuss possible extensions that 
may be important for a theoretical understanding of nuclei far from the 
stability line. 

1.2 Brief review of the standard theory 

Migdal j2] has applied Landau's theory of interacting Fermi systems to atomic 
nuclei. Here one has first to deal with two kinds of fermions — protons and 
neutrons — and second with a relatively small number of particles. In the 
past 30 years Migdal's theory has been apphed successfully by many groups 
to various nuclear structure problems. (See, for example, the reviews [3 
Ca 13 Ca Cni im and the book m.) 

Landau's theory deals with infinite systems of interacting fermions, such 
as liquid ^He or nuclear matter. If there were no interaction, the system 
would simply be a collection of independent particles, each characterized 
by its spin and a wave number k. Landau's basic assumption is that the 
interacting system can be obtained from the non-interacting one by an adi- 
abatic switching on of the interaction. In particular, there should be a one- 
to-one correspondence between the single-particle states of non-interacting 
systems and the so-called quasiparticle states in interacting systems. These 
quasiparticles behave in the correlated system like real particles in a non- 
interacting system. They obey Fermi-Dirac statistics and occupy, like the 
non-interacting particles, corresponding quasiparticle states up to the Fermi 
energy. In order to define such quasiparticles Landau considers the total 
energy of an interacting system as a functional of the occupation function 
n(k) of the quasiparticle states k. The quasiparticle energies are then given 
as first functional derivatives of the total energy with respect to n(k) and 
the interaction between the quasiparticles is defined as its second functional 
derivative. Using this approach, one can calculate properties of the excited 
system of a Fermi liquid, such as the zero sound mode in '^He. In infinite 
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systems a quasiparticle differs from a real particle essentially by its mass, 
because both can be described by a plane wave. The effective mass of the 
quasiparticle is, in general, momentum-dependent and is deduced from ex- 
periments. The renormalized quasiparticle interaction depends on spin and 
momenta and is expanded at the Fermi surface in terms of Legendre poly- 
nomials with free parameters — the well-known Landau parameters — that are 
also determined from experiments. 

Migdal extended these ideas to finite Fermi systems and applied his TFFS 
to atomic nuclei. Here the quasiparticles are the single-particle states of the 
nuclear shell model, which can be obtained experimentally from the neigh- 
boring odd-mass nuclei of closed shell nuclei. The quasiparticle interaction 
here is defined in the same way but it is isospin dependent. As in the infinite 
system, one expands the interaction at the Fermi surface in terms of Legendre 
polynomials and the parameters of the expansion, the famous Landau-Migdal 
parameters, are considered universal. They are also determined from exper- 
iment. 

The TFFS is for that reason a semi-phenomenological microscopic the- 
ory. It is a microscopic theory because all its fundamental equations are 
derived rigorously from first principles, however it also contains phenomeno- 
logical aspects such as the above-mentioned quasiparticle energies and the 
quasiparticle interaction. All the parameters are well-defined microscopically 
and in principle could be calculated starting from the bare nucleon-nucleon 
interaction. Such calculations, however, are very involved and in actual cal- 
culations one has to make severe approximations so that we can not expect 
to obtain full quantitative agreement with the phenomenological parameters. 
Nevertheless the calculated interaction parameters are in surprisingly good 
agreement with the phenomenological ones [T7] . 

1.2.1 The Landau-Migdal interaction 

As mentioned above, the energy E of an interacting system may be considered 
to be a functional of the occupation functions n(k) of the quasiparticles 
E = ii^(n(k)). If one excites the system, one basically changes the occupation 
functions by an amount (5n(k). The corresponding change of the energy is 

6E = Ek e°(k)5r2(k) + Ek k' /(k, k')6n{k)6n{k') = 

(1.1) 

Ek[e°(k) + Ek' /(k, k')6nik')]6n{k) = Eke(k)5n(k), 
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where the e''(k) are the equihbrium energies of the quasiparticles. 

The quasiparticle energies e(k) and the interaction between a quasipar- 
ticle and a quasi- hole /(k, k') defined in this way are the first and second 
derivatives of the energy functional with respect to the occupation functions: 

In nuclear matter we have an explicit spin and isospin dependence: 

/(k, k') = F(k, k') + F'(k, k')T ■ t' + [^(k, k') + G"(k, k')T ■ rV ■ a'. (1.3) 

In addition to the central interaction one has in principle also to consider 
tensor and spin-orbit forces which, however, have been neglected in most of 
the actual calculations. If one denotes the momenta before the collision by 
ki = k and \^2 = k' then, by translational invariance, the momenta after 
the collision are given by ks = k + q and k4 = k' — q, where q is the 
momentum transfer. In the third article of ref. P Landau showed that the 
main contribution (singularity) to the full scattering amplitude should come 
from small q. He also renormalized the integral equation for this amplitude 
in such a way that it contained a microscopic analog of the function /(k, k') 
and an integration over a small region near the Fermi surface. In this case 
the interaction depends only on the angle 6 between k and k'. This suggests 
an expansion in Legendre polynomials: 

F{k,k') = Y,FiPi{cos9). (1.4) 

The constants Fi are the famous Landau-Migdal parameters. One intro- 
duces dimensionless parameters by defining 

Fi = Cofi, (1.5) 

where Co is the inverse density of states at the Fermi surface: 

Co = {fr\=..- (1-6) 

Likewise, one may also expand the other components of the interaction (1.3) 
and define the parameters //, gi, g[ for the various terms. 
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Some of the Landau parameters can be related to bulk properties of the 
nucleus, such as the compression modulus 



2 d\E/A) _h^kl 



(1 + 2/0) 



(1.7) 



dk'p 2m* 



the symmetry energy /3, 



S 2m* 



(1.8) 



and the effective mass m* 



m*/m = 1 + -/i. 



(1.9) 



The various parameters have to be determined from experiments and 
can then be used to predict other experimental quantities. In particular, 
these parameters enter into the equation for 6n(k,u), the solution of which 
determines small amplitude excitation of Bose type in Fermi systems. In 
nuclei these excitations correspond to the giant multipole resonances that we 
are going to investigate. 

If one restricts the expansion in Eq. (1.4) to the lowest order, / = 0, then 
the interaction Eq. (1.3) is a constant in momentum space, which corresponds 
in r-space to a delta function in (r — r'). The next order, / = 1, is the 
derivative of a delta function, which introduces a momentum dependence 
into the particle-hole (ph) interaction. In the application of TFFS to nuclei 
nearly all calculations have been performed with only the / = component 
of the interaction, which corresponds in r-space to the following form of the 
Landau-Migdal interaction: 



F(r, r') = Co[/(r) + /'(r)Ti ■ + {g + g'r^ ■ t^)(t^ ■ cr2]5(r - r'), (1.10) 



In finite nuclei, one has to introduce density- dependent parameters be- 
cause it is obvious that the interaction inside a nucleus is different from the 
interaction in the outer region of the nucleus. The most often used ansatz is 



and similarly for the other parameters. Here po{r) is the density distribution 
in the ground state of the nucleus under consideration and fin and fex are 
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the parameters inside and outside of the nucleus. In actual calculations it 
turned out that g and g' depend only weakly on the density, so that one 
uses the same parameters inside and outside. This density dependence of the 
Landau-Migdal interaction is the basic reason for its success and universal 
applicability. 

1.2.2 General theoretical background 

Most of the results discussed so far were obtained within the phenomeno- 
logical version of the theory of Fermi liquids. In the third article of ref. P] 
Landau gave a microscopic justification of his original phenomeno logical the- 
ory using the GF technique. There he derived in a fully microscopic way the 
basic equations of his theory. He showed, in particular, that the scattering 
amplitude of two quasiparticles is connected with the response function R 
(and two-particle GF) in the ph channel and that the quantity /(k, k') is 
connected with the forward scattering amplitude of two quasiparticles at the 
Fermi surface. 

The starting point of the TFFS are the equations for the energies and 
transition amplitudes of excited states in even-even nuclei and the equations 
for moments and transitions in odd mass nuclei. As in Landau's theory, these 
equations can be also derived using the GF technique. 

The one-particle GF is defined by 

Gi2{tiM) = -t{A0\Tai{ti)at{t2)\A0) (L12) 

and the two-particle GF by 

Ki234iti,t2;h,U) = {A0\Tai{h)a2it2)atit3)atiU)\A0). (L13) 

Here a and denote time-dependent annihilation and creation operators, 
\A0) is the exact eigenfunction of the ground state of an A-particle system 
and T is the time-ordering operator. In the following we use the single- 
particle states (f\{r, s) of a nuclear shell model as basis states; therefore the 
subscripts 1...4 stand for a set of single-particle quantum numbers. 
We also need the response function R, which is defined as 

i?(12, 34) = K{23, 14) - G(21)G'(34), (L14) 

where the indices denote space and time coordinates. 
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The calculation of the excitation spectra of even-even nuclei and their 
transition probabilities can often be reduced to the calculation of the strength 
function, which describes the distribution of the transition strength in a 
nucleus induced by an external field V^{uj) 

Siu) = ^ \{An\V'^\A0)\^6{u;-u;r,), (1.15) 

where Un = — Eq is the excitation energy, while \An) and \AQ) refer to 
the excited and ground states of a nucleus with mass number A, respectively. 
As the response function Riui) in the energy representation has the spectral 
expansion 



/°° de ( 
— i?i2,34(t^,e) = XI 
-oo lT\l „ \ 



AnO* J,nO J,nOj,nO* 
^21 V^43 V^12V^34 



(1.16) 
where 

0^0 = (An|a+a2|A0) (1.17) 

is the transition amplitude between the ground state and the excited state n, 
the strength function S{uS) is completely determined by the response function 

S{uj) = - lim Im J2 <^i2,34(^ + i^)Vg. (1.18) 

^ ^^+° 1234 

This response function is defined by the Bethe-Salpeter equation in the 
ph channel [2]: 

^i2,34(^^, e) = -^31(6 + uj)G2i{e) + 

/CO rJJ 
— f/56,78(c^,e,e')^78,34(u;,eO, (1-19) 

which is shown graphically in Fig. 1.1. 

In Eq. (1.19) G is the exact one-particle GF, and U is the irreducible 
amplitude in the ph channel, which is the unrenormalized ph interaction. 
(For details see ref. [IH1-) The one-particle GF G and the response R are 
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related by a system of nonlinear equations (21 El- In particular, this system 
of equations includes the relation 

where S is the so-called mass operator related to the one-particle GF by the 
Dyson equation 

Guie) = + EG'?3(e)S34(6)G42(6) (1.21) 

34 

and G'j'2(e) = ((e — P^/2m)~^)i2 is the one-particle GF of a free particle. 



1.2.3 Microscopic derivation of the basic equations of the TFFS 

Here we will briefly describe the derivation of the basic equations of the 
standard TFFS in such a way as to clarify the explicit inclusion of complex 
configurations to be discussed in section 2. For a more detailed derivation 
see refs. 0113. 

As mentioned before, the GF G and the response function R are de- 
termined self-consistently by a system of non-linear equations. This is in 
principle an exact formulation of the (non-relativistic) A-particle problem, 
but is of little use for practical applications. 

In order to arrive at solvable equations, one applies Landau's quasipar- 
ticle concept and his renormalization procedure, which he developed for the 
microscopic theory of Fermi liquids p. For the nuclear many-body problem, 
it was done by Migdal j2] in his TFFS. 

As a first step, one splits the one-particle GF into a quasiparticle pole 
part which is diagonal in the shell model basis and a remainder, 

G,2{e) = S,2a,{^—^ + + G\,{e). (1.22) 

e — ei + to e — ei — to 

Here ai denotes the single-particle strength of the shell model pole, ei is the 
single-particle energy, n is the quasiparticle occupation number (1 or 0) and 
G^{e) is that part of the exact one-particle GF that remains when the shell 
model pole has been removed. It is assumed in Migdal's TFFS that G^{e) 
is a smooth function of e in the vicinity of ep. The dominance of the first 
term in Eq. (1.22) with respect to the e-dependence of the one-particle GF 
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is one of the basic assumptions on which the TFFS rests. The vahdity of this 
assumption is crucial for the rehabihty of the results of all TFFS calculations. 
As we will see in the following, there might be additional pole terms in the 
expansion of Eq. (1.22) that originate from more complex configurations and 
which have to be considered in addition to the simple shell model poles. This 
extension of the theory will be the main topic of this review. 

The main goal in the following is to obtain an equation for the response 
function that can be solved in practice because R{u!) contains all the informa- 
tion we need. In Eq. (1.19) the product of two one-particle Green functions 
enters. Here it is important to realize that with the quasiparticle ansatz in 
Eq. (1.22) this product can be separated at low transferred energies into a 
part A, which depends strongly on the energy u and has a 5-function maxi- 
mum with respect to e for ei approaching e2 (that is, near e^) and a weakly 
energy-dependent part B, 

- G{e + uj)G{e) = A{e, uj) + B{e, uj). (1.23) 

A is given by 

^i2,34(e, i^) = 2maia26n62A—^ — ^^(5(e - ^^-^-^) (1.24) 

ei- €2- UJ 2 

whereas B, which contains all the rest, does not give rise to a pronounced 
^-dependence. 

We now insert Eq. (1.24) into the equation for the response function and 
integrate over e to obtain in compact notation 

R{uj) = {A + B) - {A + B)UR{uj), (1.25) 

where A is the shell-model ph-propagator: 

f°° de ^ ^ ^2 — ni 
-4l2,34(c^) = / — -4i2,34(e,t^) = aia2d)i3d)24 • (1-26) 

J-oo ZTTI El — €2 — uj 

With the help of Landau's renormalization procedure one can rewrite 
Eq. (1.25) in such a way that only the known A appears explicitly in the 
equation, whereas the unknown B changes U into the renormalized ph in- 
teraction F and gives rise to effective charges for the external fields. We 
introduce a renormalized response function R{uj), which is connected with 
the original response function R{uj) by the relation 

R = e^Reg + Be^, (1.27) 
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where 

eg = l- FB (1.28) 
and the renormahzed ph interaction F satisfies the integral equation 

F = U - UBF. (1.29) 

A detailed investigation of these equations shows that F depends smoothly 
on the energy variables. For that reason the energy dependence is neglected 
when F is parameterized. The equation for R{uj) is given by 

R = A- AFR (1.30) 

and the explicit form is 

6 ( \ ^2 - ni n2 - ni ^ ~ 

J^12,MV^) — C13C24 2^ -T 12,56-^56,34- U-'J-'-J 

ei-e2-u) ei - ea - u; ^ 

This is the basic equation of the TFFS. 

Due to the renormalization procedure only the Iplh configurations appear 
explicitly in the final equation, whereas all the more complex configurations 
give rise to renormalized quantities: the ph-interaction F and the effective 
charges Cq. These quantities are not calculated within the theoretical frame- 
work but are parameterized in the way that was discussed before for the 
ph-interaction. Due to conservation laws the electric proton and neutron op- 
erators have = 1, = 0, respectively, and only spin-dependent operators 
have to be renormalized. In all practical cases the second term on the right 
side of Eq. (1.27) does not contribute, so it is sufficient to solve Eq. (1.31). 
Here all quantities are known. The single-particle energies are given by the 
nuclear shell model or taken, as far as possible, from experiment and the 
parameters of the interaction and the effective operators have been deduced 
from experimental data. 

Equation (1.31) has exactly the same form as the conventional RPA equa- 
tion that, however, has been derived using approximations from the outset. 
In the present derivation, Eqs. (1.27) and (1.28) are still exact. For that 
reason one is able to obtain relations (the Ward identities) between the ef- 
fective operators and the effective interaction. In cases where conservation 
laws exist these relations determine the effective operators completely. In 
addition, the present derivation shows more clearly the range of validity of 
the theory, which naturally also applies to the conventional RPA equation. 
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Moreover, the GF formalism provides a natural basis for an extension of the 
theory, as we will see in the next section. 

For practical reasons one solves not Eq. (1.31), but the related equation 
for the change of the density matrix puiw) in the external field V^{u), which 
is defined as 

Puiuj) = -J2Ru,34i^)egVgiu;). (1.32) 

34 

The equation for pi2 follows from Eq. (1-31) and has the form 

Pl2{uj) = -J2 ^12,3465 V;°3 - J2 ^12,34^34,56P56(^^)- (1-33) 
34 3456 

The expression for the strength function is then given by 

S{uj,A) = -- Im ^e,Oi2(^ + ^A), (1.34) 

^ 12 

where A is a (finite) smearing parameter which simulates a finite experimen- 
tal resolution and at the same time phenomenologically can include configu- 
rations not dealt with explicitly in the approach under consideration. 

Equations (1.33) and (1.34) are the main equations that are used in the 
calculations within the TFFS. 

1.2.4 Coordinate representation 

Most of the calculation that we will present here have been performed in 
r-space and not in the configuration space of a shell model basis. For the 
cases of RPA and TFFS the method was suggested in refs. [201 iU] and was 
included by us in the ETFFS [TS]. The main reason for this choice is that 
the r-space representation is much more appropriate for the treatment of the 
single-particle continuum, as first pointed out in refs. [201 1^- Therefore we 
give here some relevant equations in the coordinate representation. 
Equation (1.33) has the form 

p(r, ^) = - / ^(r, r', cj)eg1^°(r', uj)d^r' - (1.35) 
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The ph propagator A, given by 



A{r,r\uj) = j: ^l{r)Mr'W*2{rWi{r% (1.36) 

^ ei - e2 - u; 

can be rewritten as 

A{r, r', u;) = -Y^ m^l{r)y^i{r') [G{r', r; + u) + G{r', r; e^ - u)] (1.37) 
1 

using the formula for the one-particle GF 

G(r.r';.) = i:^2W^. (1.38) 

2 e - £2 

where v?2(i") are the single-particle wave functions calculated in a mean-field 
potential. 

The summation in Eq. ()1.37p is over states below the Fermi surface, i.e. 
the single-particle continuum is already contained in Eq. (1.38). On the 
other hand, the coordinate part of this GF Gij = gij/rr' can be expressed 
in closed form in terms of the regular Pij^ and irregular ylj^ solutions of the 
one-dimensional Schrodinger equation as 

9iAry;e) = |?l/;;^(r<;6)|/(f (r->;6)M,(6), (1.39) 

where r< and r> denote the lesser and the greater of r and r', respectively 
and W is the Wronskian of the two solutions. The irregular solution i/ij^ is 
determined by the boundary conditions at oo; e.g., for neutrons 

/ox 

Uij \^ ~^ cxd) ~ exp{—kr) (1-40) 
for negative energies e < and 

y^if{r ^ cx)) ~ exp[i{kr - y + ^ij)] (1-41) 



for positive energies e > 0, where k = ^J2m\e\/?i and 5ij is the scattering 
phase for the mean nuclear potential considered. 

Thus, the functions ' are calculated numerically if the mean potential 
is known. For e < the functions gij have no imaginary part; that is, the 
Iplh states have automatically no width if the smearing parameter A = 0. 
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Inclusion of the single-particle continuum makes it possible to obtain a 
physical envelope of the resonance without using a smearing parameter, that 
is, to obtain directly the escape width F^. 

In the complex configuration problem under consideration, using the rep- 
resentation of the single-particle wave functions (A-representation) gives ma- 
trices of a very large rank especially, in the case of treatment of the ground 
state correlations caused by complex configurations. Using the coordinate 
representation affords a big numerical advantage in this problem because the 
rank of matrices is determined not by the number of configurations but by the 
number of mesh points used in solving the corresponding integral equation. 

1.2.5 Main results of the GMR description within the TFFS 

In order to distinguish GMR from other collective excitations, one can define 
them as follows: 

1. The form and the width F of the resonance depend rather weakly on 
A; as a rule, the dependence F ~ is used. 

2. The resonance mean energy E also depends weakly on A; usually one 
uses E ~ 

3. The resonance width is small compared with its excitation energy. 

4. The resonance exhausts a large fraction of its energy weighted sum rule 
(EWSR) — usually more than 50 percent. 

The last of these is the most quantitative characteristic of the resonance and 
justifies its name "giant" . 

In 1971-72 in inelastic electron j23 122] and proton [21] scattering, giant 
multipole resonances (GMR) were detected that were different from the well- 
known isovector electric dipole resonance. That was the starting point of a 
period, sometimes called a renaissance of giant resonance physics, of very 
rapid and intensive development. A large amount of experimental data on 
the GMR in stable nuclei — principally their energies, total strengths, widths 
and resonance envelopes — has been accumulated and discussed since the mid- 
seventies. Currently there exists experimental information on more than 20 
different types of GMR that were detected in a large number of nuclei in 
a broad range of excitation energies. A detailed review, with experimental 
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results up to the end of the eighties, is presented in ref. j2Sl- Further informa- 
tion can be found in the proceedings of the last four international conferences 
devoted to giant resonances [211 I2H1 1211 , and in a recent monograph jSU] 
that gives an excellent review of the present experimental situation. Conven- 
tional theoretical methods such as the RPA and QRPA, and their comparison 
with the data, have been the subject of many review articles and books (e.g., 
refs. jSniinSll2SllH3 and [HZl EOl 1121 El EHl El IE] ) • For nuclei with A > 40, 
the experimental situation is essentially settled. There is no longer any ma- 
jor controversy over the centroid energies and the total strengths of GMR, 
and the theoretical interpretation of these data within the framework of the 
conventional theoretical methods is also clear. These methods, however, do 
not allow description of the widths and the fine structure of the resonances, 
nor do they offer any possibility to analyze complex spectra with overlapping 
resonances. This is the main subject of our review and will be presented in 
the following sections. 

In the following we illustrate and briefly discuss typical results that have 
been obtained within the continuum TFFS (CTFFS), which is, as men- 
tioned before, formally identical to the continuum RPA (CRPA). We show 
in Figs. 1.2-1.4 and Table 1.1 the hadron strength functions for the E2 
isoscalar (IS) and isovector (IV) resonances in ^°^Pb and the E2 IS reso- 
nance in ^°Ca and some of their integral characteristics. The calculations 
have been performed in the coordinate representation within the CTFFS, 
i.e., using Eqs. (1.35), (1.37) and (1.39). In order to simulate the finite ex- 
perimental resolution, we introduced a smearing parameter A with a value of 
250 keV. As the TFFF is not self-consistent, as discuss in the beginning, we 
have to determine the parameters of the effective interaction, Eqs. (1.10) and 
(1.11), from experiment by fitting some specific theoretical results to exper- 
imental data. In our calculations we always used, with the exception of f^x, 
the following Landau-Migdal interaction parameters, which were adjusted 
previously to various experimental quantities |10H I102j : 



fin 


= -0.002 


f 

J ex 


= 2.30 


f 


= 0.76 


9 


= 0.05 


9' 


= 0.96 


Co 


= 300 MeVfm3 



(1.42) 
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For the nuclear density po(^) in the interpolation formula (1.11) we chose 
the theoretical ground state density distribution of the corresponding nucleus, 

PoW= E ^(2ji + lR'(r), (1.43) 

where Ri{r) are the single-particle radial wave functions of the particular 
Woods-Saxon potential used. For other details of the calculations, including 
the definitions used in Table 1.1, see section 3.1.3 

For the parameter /ex we have used the values fex = —1.9 and -2.2 for 
^°^Pb and ^°Ca, respectively. These parameters were adjusted to reproduce 
the energies of the first excited 2^ level in 208p|^ ^^-^^ g- [gyg[ ^''Ca. We 
will see in sections 3 and 4 that the same parameters (jl.42p can also be used 
if more complex configurations are considered explicitly. Even the parameter 
fex changes only slightly in the latter case. We have used these parameters 
in all our calculations for stable and unstable closed shell nuclei from ^^O 
to ^°^Pb, where we investigated many different types of GMR. In our opin- 
ion, the reasonable agreement with experiment we obtained confirms the 
assumed universality j2] of the parameters of the Landau-Migdal interaction, 
Eq. (1.10). 

In Figs. 1.2-1.4 the CTFFS results are given together with the results of 
the calculations without taking the effective interaction into account ("free 
response"). As the parameters fin, fex ifLyfL) isoscalar (isovector) 

resonance are negative (positive), the IS (IV) resonances are shifted to lower 
(higher) energies, when the interaction is included, compared with the free 
response. For the latter we have the shell model estimate for 208p|-| ^ ^ 
2 X 41A^i/3 = 13.8 MeV while, according to Table 1.1, Ei^ = 8.1 MeV, = 
18.2 MeV or, if one uses another definition for E = i?2,o (see section 3.1.3), Eis 
= 10.2 MeV, Ei^ = 19.3 MeV. In both cases — with and without interaction — 
reasonable depletion values (90-100 percent for large energy intervals) of the 
corresponding EWSR have been obtained (see Table 1.1). The depletion 
is in satisfactory agreement with the corresponding experimental values of 
the EWSR. (See sections 3 and 4 and, in particular. Tables 3.3 and 3.7.) 
In Table 1.1 we give also the experimental values of the mean energies. It 
should be noted, however, that the experimental data were obtained, as a 
rule, for intervals which are much smaller than those given in Table 1.1, not to 
mention that the experimental mean energies may be determined in different 
ways that may be important for unstable nuclei for which taking the single- 
particle continuum into account is necessary. Therefore the comparison with 
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experiment should be made more carefully, and this will be done in section 3, 
where we will also discuss the width question. 

Thus, we have obtained a relatively good description of the mean energies 
and total strengths of the GMR under consideration. These are typical results 
of the CRPA. 

In Figs. 1.2-1.4 we have chosen the smearing parameter noticeably larger 
than the experimental resolution in order to simulate at least some of the 
decay width not included in the present approach. Nevertheless, one cannot 
see in Figs. 1.2-1.4 any resemblance to an observed resonance because the 
smallest width among the three resonances under discussion is the one for the 
E2 IS resonance in ^°^Pb, which has an experimental width of 3.1 ± 0.3 MeV 
j35j . In other words, the widths of resonances, which are among the most 
important characteristics, are not reproduced within the CRPA. Even in 
medium-mass nuclei, where the role of the continuum (escape width)is much 
larger, the theoretical widths are still in disagreement with the experimental 
ones. The reason is well known: the one particle-one hole configurations 
describe only the escape widths, which — in general — are only a small frac- 
tion of the total widths. For a realistic description one has to include the 
spreading widths, which are due to the coupling of the one particle-one hole 
configurations to more complicated configurations. One possible solution to 
this problem will be discussed next. 

1.3 Physical arguments for extending the standard ap- 
proach 

We have seen that the Landau-Migdal theory is based an a microscopic many- 
body theory with, however, important elements taken from experiment. For 
that reason it is quite natural that an extension of that successful theoretical 
frame work is also based on experimental facts. It is well known from nuclear 
spectroscopy that in odd-mass nuclei that differ from a magic nucleus by 
one nucleon (or hole), the coupling of the low-lying phonons of the even 
nucleus to the odd particle (or hole) play an important role. It gives rise to 
a strong fragmentation of the corresponding single-particle (hole) strength 
over a range of the phonon energy. This observation gives us the possibility 
to include in the standard TFFS the coupling of the Iplh states to more 
complex configurations. We will show that with this coupling to the low-lying 
phonons one includes those 2p2h configurations that give rise to the strongest 
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fragmentation of the giant resonances, which — together with the couphng to 
the single-particle continuum — makes it possible to calculate quantitatively 
the strength distribution of the GMR. 

1.3.1 Giant resonances in nuclei: the width problem 

It is clear now that GMR are a universal property of nuclei. The investi- 
gations of GMR are not only important for a detailed understanding of the 
structure of nuclei, but they are also an important tool for a better under- 
standing of nuclear reaction mechanisms involved in the excitation of the 
different types of these resonances. Moreover, we obtain from the investiga- 
tion of GMR additional information on the Landau parameters. The most 
important of these is fin, which is connected with the compressibility of nu- 
clear matter and is therefore of crucial importance in astro-physics. The 
isoscalar electric monopole resonance (breathing mode), on the other hand, 
is closely related to this parameter and therefore one would need to know this 
resonance in nuclei far from stability in order to obtain the dependence of the 
compression modulus on the number of protons and neutrons. The magnetic 
resonances are related to the parameters g and g' of the spin- and isospin- 
dependent parts of the forces. The latter is related to pions in nuclei and 
is of special interest in connection with the possibility of pion condensation. 
Phenomenologically, GMR inform us about the nuclear shape (splitting of 
the El, E2and EO resonances in deformed nuclei), and about volume, surface 
and other kinds of vibrations. 

The understanding of the widths of GMR is obviously connected with 
the damping of small amplitude vibrations in finite systems, as we shall 
soon discuss. Thus the general problem of how energy from highly ordered 
excitations is dissipated in nuclei, including the question of transition from 
order to chaos, can be clarified through GMR studies ^Tj . 

Many ideas from GMR physics have been used in other applications, 
such as the recent investigations of metallic clusters [IHI and the fuUerene 
molecules jSU]. See also ref. jSHI, chapt. 11. 

Necessity of inclusion of complex configurations 
and single-particle continuum 

As discussed in section 1.2.5, the standard continuum TFFS or the con- 
tinuum RPA in closed shell nuclei are able to describe only two integral char- 
acteristics of GMR: their mean energies and total strengths. The quantity 
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that is of equal importance, the strength distribution of the GMR (i.e., their 
widths) cannot be reproduced within this approximation. The reason for this 
failure has been already been indicated above: the coupling of the Iplh con- 
figurations to more complex configurations, which gives rise to the spreading 
width, has so far been neglected. The escape width, which is included in the 
present continuum approaches, represents in the giant resonance region only 
the lesser part of the total width. There are several reasons why one would 
like do describe theoretically the widths of the GMR quantitatively: 

1. the intellectual challenge to develop a microscopic theory that gives a 
quantitative explanation for the collective motion in strongly interact- 
ing finite Fermi systems; 

2. the new insight into the fine structure of the GMR due to the rapid 
improvement of the experimental resolution to AE < 10 keV, which 
needs to be understood (see, for example, ref. 

3. the need for microscopically derived strength distributions that quanti- 
tatively reproduce the date on resonances in the medium mass region, 
where the various multipole resonances overlap in energy. The con- 
ventional analyses with phenomenological transition densities are no 
longer applicable because they introduce strong uncertainties, as we 
shall discuss in the following sections. 

The inclusion of the single-particle continuum gives the physical enve- 
lope for processes at excitation energies higher than the nucleon separation 
energy. For giant resonances this gives the escape width F^^, the magnitude 
of which depends significantly on the mass number, the excitation energy, 
the multipolarity, etc., so that it is not justified, especially in the calcula- 
tions of such delicate properties as fine structure and decay characteristics, 
to simulate the role of the continuum by a constant smearing parameter. 
The exact microscopic treatment of the continuum is therefore crucial to a 
realistic theory of giant resonances. 

In addition to this, a realistic microscopic theory of collective motion in 
nuclei has also to consider more complex configurations than those included 
in the RPA. There are new data, e.g. 1144111481 1145j . and the largely unsatis- 
factory explanation of the older results concerning the low-lying structures in 
cross sections in a wide excitation energy range around the nucleon binding 
energy [S^IMll^l^llljl l^. see also section 3.2.5. An extended theory will 
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also have implications for the interpretation of experimental data obtained 
with modern germanium detectors and gamma spectrometers such as EU- 
ROBALL cluster, EUROBALL and others [Hll^. Unprecedentedly high 
resolution and high efficiency of detecting gamma rays with energies up to 
20 MeV have already given new and very precise information, not only on 
deformed nuclei, but also on low-lying levels in odd and even-even spherical 
nuclei. In fact, these detectors give direct information about the low-lying 
complex configurations containing phonons |^ I59| IHUl IHT] I62j . which may be 
seen again in the fine structure of the giant resonances. At last, in order to 
explain the large amount of available data on the decay properties of GMR's 
gained in experiments with coincidences of secondary particles |^ IHU] it is 
also necessary to take complex configurations into account (see, for example, 
refs. [HlZni). 

It is clear that in the immediate future the number of such data will in- 
crease rapidly and that these results require improved microscopic approaches 
for their interpretation. 

Ground state correlations caused by complex configurations 

The ground state correlations (GSC) problem has a long history. (See, 
for example, references in [^31111] and also the article •) It started with 
the Hartree-Fock approximation (HP), where the effects of the Pauli princi- 
ple was included in the calculation of the ground state of fermion systems. 
Some specific ground state correlations are taken into account if one cal- 
culates excited states within the RPA. The most important consequence of 
these ground state correlations is that the energy-weighted sum rule for the 
transition probabilities is conserved. This is not so in the Tamm-Dankoff 
approximation, which starts from the uncorrelated HP ground state. Dur- 
ing the past ten years, in connection with the development of the extended 
TPPS, where configurations beyond the Iplh states are included in the ex- 
cited states, the question arose as to how far one has also to consider the 
same configurations in the ground state. |HS1 IHH Ei] There is a fundamental 
difference between the effects of GSC in the RPA and their effects in models 
where more complex configurations are included explicitly. The RPA GSC 
do not lead to the appearance of new transitions compared with the TDA, 
but only shift the energies and redistribute the transition strengths. 

The GSC induced by the more complex configurations, on the other hand, 
lead not only to a redistribution of strength but also to new transitions. 
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which give rise to a change of the EWSR jHEl ESI- Thus these GSC are 
at least as important and physically interesting as the RPA GSC. Actually, 
their consequences are much richer and far-reaching. The present approach 
is an extension of the previously developed Extended Second RPA (see, for 
example refs. |96llS^ ). in which uncorrelated 2p2h GSC have been considered. 

We will see in the following that these effects play a noticeable — sometimes 
decisive — role in the theoretical description of the experimental data. The 
most striking example obtained within the GF approach is the explanation 
IHH] of the observed Ml excitations in '^''Ca and ^^O with energies of 
about 10 MeV and 16 MeV, respectively, solely as a result of ground state 
correlations. 

1.3.2 Nuclei far from the stability line 

There is increasing interest in the structure of nuclei far from the stability 
line. The study of these exotic nuclei, is of importance not only in itself 
1111121 ED, but also for its relevance to astrophysics jH |13 IH IH EI] • The 
ETFFS we are discussing here may play an important role in the analysis of 
the experiments done at the proposed radioactive beam facilities. As men- 
tioned earlier, one of the crucial quantities one wants to know is the breathing 
mode in nuclei with very different numbers of protons and neutrons. This 
will give us the compression modulus as a function of the proton and neutron 
number, which is needed for the extrapolation to nuclear matter. We may 
suppose that in nuclei far from stability, even with closed shells, the high- 
lying spectra may be as complicated as in the medium mass nuclei, where the 
various multipole resonances overlap and a microscopic theory is necessary 
for the analysis in order to obtain reliable nuclear structure information. 

If we extrapolate our present knowledge of unstable nuclei to nuclei far 
away from the stability line we may expect two characteristic features: (i) 
there will be very low-lying collective states and (ii) the nucleon separation 
energy may also be relatively low. For these reasons a realistic theory has to 
treat the continuum in an exact way, and the phonon coupling is not only 
important for the analysis of GMR but also for a quantitative understand- 
ing of the low-lying spectrum. In connection with the application of the 
present approach, one has to investigate the extent to which the Landau- 
Migdal parameters are dependent on the numbers of protons and neutrons. 
Our extensive experience indicates that this dependence may be quite weak, 
so we can at least use the present parameters as a good starting point. The 
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second important input into the theory concerns the single-particle spectrum 
and the single-particle wave functions. Here one may use self-consistent ap- 
proaches, for example, those with the density functional, in order to obtain a 
reliable quasiparticle basis. Investigations in these directions are in progress. 

It should be emphasized that only a reliable inclusion of the single-particle 
continuum can make it possible to do calculations for nuclei with the separa- 
tion energy near zero. This is important for understanding drip- line nuclei and 
for astrophysical sdudies. For neutron-reach nuclei with the separation en- 
ergy near zero, this is of prime interest because of the absence of the Coulomb 
barrier. The CRPA calculations in ^^O have shown that the strength distri- 
butions of the E2 jin] and isovector El j37| resonances are very different from 
those for ^^O; the resonances are more spread out, shifted down and have a 
noticeable low-lying strengths. The effect of complex configurations is also 
noticeable, at least for the isoscalar E2 resonance in ^^O |1H]- It should be 
pointed out, however, that except for ref. jlH] and the calculations we shall 
present in section 3.3, there exists almost no theoretical information about 
the role of complex configurations in unstable nuclei |3Uj . 

1.3.3 Implications of experiments for and current status of the 
microscopic theory 

We can summarize our discussion so far by asserting that a microscopic 
theory that is able to describe quantitatively the structure of collective exci- 
tations in nuclei and which is based on the mean-field approximation has to 
include four major effects: 

1. the Iplh RPA, which creates collectivity out of the uncorrelated ph 

states, clS db starting point; 

2. complex configurations beyond the Iplh states, which give rise to a 
fragmentation of the collective states derived from the RPA; 

3. the single-particle continuum; 

4. ground state correlations induced by the complex configurations under 
consideration. 

In addition, one should not use separable forces, because one then must use 
different forces for each mulipolarity, which strongly reduces the predictive 
power of the theory. Indeed one needs an interaction that is universal for 
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the whole periodic table, or at least that changes only very little with the 
mass number, and which should be adjusted to quantities other than those 
that one is going to calculate. As we shall see, the GF approach that we are 
going to discuss in the what follows allows the inclusion of all these effects 
simultaneously. 

It is obvious that, compared with the simple Iplh configuration problem, 
the present task is much more difficult — both theoretically and numerically. 
In addition we shall develop and apply various stages of sophistication of 
our theory to the nuclear structure problem in order to clarify the different 
effects. At present there exist several other approaches that have considered 
some of the effects mentioned above. These are reviewed in refs. |4H 15^ 
("pure" 2p2h configurations) and jHSHSnHnO] (configurations with phonons). 

In the past, microscopic theories of GMR have been developed using two 
different approaches: RPA + continuum on the one hand and RPA + complex 
configurations on the other. As we have seen, however, both extensions of the 
RPA are need to explain the data. The first of these can now be considered 
solved, and there exist several numerical methods for it. One, which was 
mentioned in section 1.2.4, uses the GF method. There one considers the one- 
particle continuum exactly (for a contact interaction) by transforming the 
RPA equation into the coordinate representation. Other methods for solving 
this problem have also been developed that even admit the use of nonlocal 
forces jTU [TO] . As for the problem of including complex configurations, the 
most advanced approach is the quasiparticle-phonon model for magic and 
non- magic nuclei by Soloviev and his co-workers jHH EHl- These authors, 
however, used separable forces in order to reduce the numerical difficulties of 
the problem, and they leave out the single-particle continuum. In addition, 
the ground state correlations are included only partially, that is, mainly on 
the RPA or QRPA level. 

The microscopic theory for GMR that satisfies the requirements men- 
tioned above turned out to be quite difficult to formulate and especially to 
realize numerically if one uses non-separable forces, as shall will do. 

There have been some successful developments in this direction in the past 
ten years. The first attempts of this kind, which simultaneously consider RPA 
configurations, the single-particle continuum (escape width F |) and more 
complex (2p2h [721 [73] or lplh®phonon [721 [Hj) configurations (spreading 
width F I) using non-separable forces were made in refs. [721 [731 [721 [H] for 
some closed shell nuclei. These authors investigated various types of GMR 
using, of course, different approximations and methods. 
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The model developed in ref. [73] used only lplh®plionon configurations 
and it considered only a particle-phonon interaction. It is based on the Bohr- 
Mottelson model for the strength function of the phonons. This model was 
also successfully used [7^] to calculate partial branching coefficients of the 
proton decay of the isobar-analog and Gamow-Teller resonances in ^°^Bi. 
The papers [73J [70] and [72] were the first articles in which all three (that 
is, the above-mentioned items 1, 2 and 3) microscopic mechanisms of GMR 
formation were used to explain such delicate properties as the decay char- 
acteristics of the GMR. It was also shown in refs. [7^ [TH] that the complex 
lplh®phonon configurations noticeably improve the description of the decay 
characteristics. The advantage of this method [731 IZH] is a self-consistency 
(on the RPA level), that is, the phonons that are used in the extended theory 
have been obtained in RPA using the calculated interaction. In this develop- 
ment, however, ground state correlations due to complex configurations have 
been ignored. 

The most extensive investigations of GMR, which include the above- 
mentioned effects, were performed within our ETFFS approach, where cal- 
culations for stable and unstable closed shell nuclei have been made. The 
theory is based on the consistent use of the GF method [73 [73 [HOI EH E2 ESI 
IH31 IHSl IHH1IHI] • The ETFFS simultaneously takes into account Iplh, complex 
lplh®phonon configurations, the single-particle continuum and ground state 
correlations both of the RPA type and of those caused by the complex con- 
figurations under consideration. In addition, in its final equations it includes 
explicitly both the effective particle-hole interaction and the quasiparticle- 
phonon interaction. 

1.4 Qualitative discussion of the extension 

A theoretical approach that takes into account the 2p2h configuration space 
including the full 2p2h interaction is numerically hardly solvable if one also 
uses a realistically large configuration space. For that reason the main ap- 
proximation in ETFFS concerns the selection of the 2p2h configurations. In 
our approach, guided by experimental observations, we include the most im- 
portant correlations in the 2p2h space by coupling phonons (correlated Iplh 
states) to a one-particle and one-hole state. With this procedure we obtain 
effectively lp0phonon and lh®phonon configurations. If we then couple an 
additional hole and particle, respectively, to the previous configurations we 
obtain lplh®phonon configurations where part of the 2p2h interaction is 
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included. As one can see from the applications, these configurations are in- 
deed the most important ones for the understanding of the spreading width 
of GMR. The configurations with phonons also nicely explain a part of the 
low-lying spectrum in the neighboring odd mass nuclei, configurations with 
phonons are used in many theoretical approaches [SHI ESI EH] ) • 

There is, however, an additional fact that greatly simplifies the problem, 
and that is the existence of a small parameter for closed shell nuclei |88j : 



where (l||g||2) is a reduced matrix element of the amplitude for low-lying 
phonon creation with the energy ujg and 1 represents the set of single-particle 
quantum numbers rii, /i and ji for spherical nuclei. Henceforth, when we 
refer to the g^ approximation, it will be understood that the dimensionless 
a is small. Using this small parameter affords following advantages: 

1. We obtain a general principle for selecting terms: as a is small, we may 
restrict ourselves to lplh®phonon configurations, which correspond to 
second order in g (two-phonon configurations correspond to terms of 
order g^). Because we use the g^ approximation in the propagators 
of our integral equations, our approach is not the usual perturbative 
theory in g^. 

2. For the widths of GMR the most important contribution comes from 
low-lying phonons, which give rise to a strong energy dependence in the 
energy range of the high lying collective (Iplh) RPA solutions. There- 
fore we may confine ourselves to the most collective low-lying phonons, 
which are restricted in number. The effects of the other phonons are 
effectively already included through the phenomenological parameters 
of our approach. 

3. The restriction to a small number of collective phonons noticeably re- 
duces the numerical difficulties. This is especially important for the 
present approach, in which non-separable forces and the GSC induced 
by complex configurations are considered. 

4. As some of the lplh(8)phonon configurations are treated explicitly in the 
ETFFS, one should expect that the Landau-Migdal parameters that are 
determined within the Iplh approximation may change. As we restrict 



(l||g||2)^ 



< 1 




a = 
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ourselves to the approximation and the collective low-lying phonons, 
this effect in the actual calculations is found to be small. 

The ETFFS approach is, like the original TFFS, a semi- microscopic the- 
ory. As our approach is based on the TFFS, we actually do not need addi- 
tional experimental input beyond that used already in the TFFS. We must, 
however, "correct" some of the parameters in order to avoid double counting. 
Such corrections can be performed fully consistently within our approach. 
The most important corrections refer to the single-particle energies which 
are taken — as far as possible — from experiment, or else from a shell model 
potential that is carefully adjusted to the corresponding closed shell nucleus. 
The single-particle wave functions are also taken from that model. These 
quantities contain contributions from the same phonons that enter the com- 
plex configurations under consideration. In order to avoid double counting 
due to these phonons, the single-particle model has to be "refined" from 
this mixing. The procedure for this will be described in section 3.1.2. The 
complex configurations that we treat explicitly in our extended approach are 
also included implicitly in the force parameters of the original TFFS ap- 
proach. Therefore, in principle, one has also to correct the Landau-Migdal 
parameters. 

There exists so far no self-consistent theoretical approach that includes 
all the effects discussed above. In such an approach one would start with 
an effective two-body interaction that would allow to determine the single- 
particle energies and wave functions and the ph interaction. As within such a 
procedure the phonon effects are not included, our extended theory would be 
the natural formalism in which to do it. In our extended version of the TFFS 
we did not include the so-called tadpole graphs with the low-lying phonons 
under consideration that have been used in the self-consistent version of the 
TFFS IB]. Their contribution is contained effectively in our "refined" mean 
field. 
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2 Framework of the extended theory 



2.1 General description of collective excitations, 
including the particle-hole and quasiparticle- 
phonon interaction 

As mentioned before, the original TFFS allows to calculate only the cen- 
troid energies and total transition strength of giant resonances because the 
approach is restricted to Iplh configurations. In order to describe more de- 
tailed nuclear structure properties one has to include higher configurations. 
Here we describe the derivation of the main ETFFS equations. These equa- 
tions contain both the quasiparticle-phonon interaction and the effective ph 
interaction in a general form. 

2.1.1 General relations 

Equation ()1.19p can be considered a definition of the response function R 
only if the quantities S and U are known. In order to obtain realistic numer- 
ical results, however, we have at our disposal only model approximations of 
these quantities because it is not possible to solve the whole system of non- 
linear equations for i?, S and U for a realistic case |2l El • For that reason 
Landau introduced phenomenological elements into his microscopic theory. 
In Landau's approach S is irreducible in the one-particle (one-hole) channel 
and U is irreducible in the ph channel. As the strongest energy dependence 
in Eqs. (1.19) and (1.21) is considered explicitly through the one-particle 
(one-hole) and the ph propagators, respectively, the irreducible parts are 
weakly energy dependent and are parameterized in an energy-independent 
way. In our extended theory the complex configurations give rise to a strong 
energy dependence in the previously weakly energy-dependent irreducible 
quantities S and f/, which we now have to consider explicitly. In order to 
do so we represent S and f/ as a sum of two terms, in which the first terms 
are again assumed to depend only weakly on the energy. As we shall see 
in the next section, these terms are irreducible with respect to the complex 
configurations that are considered explicitly and can, for the same reasons 
as in Landaus original theory, be parameterized. 




(2.1) 
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Using Eq. (2.2) one can transform our main equation, Eq. (|1.19|) . into the 
symbolic form 

R = R^- R^UR, (2.3) 
where the quantity R'^ satisfies the equation 

= -GG + GGU^R\ (2.4) 

Equation (2.3) is more convenient in the sense that it contains only the 
energy-independent amplitude f/, which allows us to apply the known renor- 
malization procedure j2] to it. 

Furthermore, we rewrite the Dyson Eq. (1.21) as 

Gi2(e) = Gi2(e) + E G13 (6)2^4(^)^42(6), (2.5) 
34 

where G is the solution of the Dyson equation with the mass operator S, 

Gi2(e) = G°2(e) + EG°3(e)E34G42(e). (2.6) 

34 

The quantities denoted with tildes, S and U , as mentioned before, are 
the analogs of the quantities S and U in the standard TFFS. The energy- 
dependent quantities Tf" and U'^ introduce the effects of the complex configu- 
rations into the formalism. The corresponding single-particle basis {(^a, ^a}, 
which is defined by E, is a new or "refined" basis and should be obtained 
from the basis {ip\, e\} used in the TFFS. 

2.1.2 Renormalization of the general equation for the response 
function and the equation for the density matrix. 

Equation (2.3), as well as Eq. (1.19), contains summation and integration 
over all states, including those that are far away from the Fermi surface and 
which, in actual calculations can be taken into account only effectively. For 
that reason it is necessary to perform the Landau-Migdal renormalization 
procedure in order to obtain an equation with the summation restricted to 
the vicinity of the Fermi energy. This is performed in analogy with the 
renormalization of the equation for the response function within the standard 
TFFS described in section 1.2.3. [TH] . 
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In analogy with Eqs. (1.25) and (1.26), we introduce the generahzed prop- 
agator A, which will explicitly be given by the complex configurations under 
consideration and the quantity B, which contains all the rest: 

/oo 
— i?^2,34(^, e) = ^12,34(^) + 5l2,34. (2.7) 
-oo zm 

It is assumed that the quantity B depends only weakly on uj compared to A. 
One obtains from Eq. (2.3) the renormalized response function R{uj), which 
is connected with our initial response function by the relation 

-Rl2,34 = X! (^g )l2,56-R56,78('^)(e5)78,34 (2.8) 
5678 56 

where R{uj) satisfies the equation 

-Rl234(t^) = ^12,34(t^) - ^12,56 (^)-^5678-R7834(t^). (2.9) 

56,34 

Here the new effective charge operator and the effective ph interaction 
amplitude F are given by 

(Cg) 12,34 = — X! -^12,56-856,34 (2.10) 

56 

-^12,34 = 51(^9)12, 56^^56,34 (2.11) 
56 

Equation ()2.4|) still contains the full GF and, in fact, according to Eq. (2.7), 
determines the desired propagator A[u}) in Eq. (2.9). It is useful to trans- 
form it to a more convenient form. Acting on both sides of Eq. ()2.4j) with 
the operator {G~^G){GG~^) and using Eq. (2.5), we obtain, in the time- 
representation, 

i?"(12, 34) = R%12, 34) + ^ i?°(12, 5Q)W^{5Q, 78)i?^(78, 34), (2.12) 

5678;t 

where 

^°(12, 34) = -^(3, 1)G(2, 4) (2.13) 
1^*^(12, 34) = PVo (12, 34) - zS'^(31)S"(24) (2.14) 
W^{12, 34) = t/"(12, 34) + zS"(3, l)G'-^(2, 4) + tG~\3, 1)2^^(2, 4). (2.15) 



31 



Because the quantity B only weakly depends on the energy as compared 
with A, it follows from Eq. (2.8) and Eq. (1.18) that S{E) = hmA_.+o S{E, A), 
with 

SiE,A) = - Im ^(e~g\/°);ii?i2,34(^ + ^A)(g,\/°)43 (2.16) 

^ 1234 

(e,ni2 = E(e.)2i,43V^34- (2.17) 

34 

We see that the strength function is completely determined by the renormal- 
ized response function R{u;). 

The advantage of these transformations is that, as in the original TFFS, 
the experimental quantity S{E) is now connected with the renormalized re- 
sponse function, which can be calculated from Eq. (2.9). In contrast to 
Eq. (2.4) for the full response function R, which cannot be used for nu- 
merical calculations, the equation for the renormalized response function, 
Eqs. (2.9), (2.16) contain only quantities which will be defined in the next 
section: the generalized propagator A, which can be explicitly calculated 
from the complex configurations, and the effective ph interaction F, and the 
effective charge operator that are treated in the same way as in the original 
TFFS approach. 

As in the case of the TFFS (section 1.2.3), we also introduce the change 
of the density matrix pi2 due to an external field for the generalized case 
under consideration: 

Puiuj) = - E^1234(^)(g,V^°)43. (2.18) 

34 

The equation for pi2 follows directly from Eq. (2.9) and has the form 

Puiuj) = Ai2,34(cu)(egV^°)43 - E ^12,34 (^)^3456P56(C^), (2.19) 

34 3456 

and the expression for the strength function can be written as: 

S{u, A) = Im Y.i^gV%,puiuj + ^A). (2.20) 

^ 12 

Equations ()2.20|) and ()2.19|) are our main general results. We stress that 
in Eqs. (2.1), (2.2), (2.4) and (2.14), (2.15) we did not yet specify the form 
of our energy-dependent terms and f/^ and the transformation from R to 
R is formally exact. So far we have only assumed that the quantities B and 
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U, as in the standard TFFS, depend only weakly on the energy. The same 
holds for iq and F. After we specify the mass operator and the irreducible 
amplitude in section 2.2, we shall obtain the generalized propagator A in 
section 2.3. With this information we are able to apply the ETFFS to GMR, 
the results of which we shall discuss in sections 3 and 4. 

Equations (2.8) and (2.9) for our renormalized response function R{uj) 
and the final equations, Eqs. (2.20) and (j2.19|) . have a structure similar to 
the relations in the TFFS given in section 1.2.3. Let us point out the main 
differences as compared with the TFFS. 

First of all, the generalized propagator A, will differ the most from the 
TFFS form, Eqs. (1.23) and (1.26), because we must now solve an equation, 
either (2.4) or (2.12) to obtain the propagator. In the standard theory one 
considers the propagation of a particle-hole pair, whereas in the extended 
theory the propagator has to include also the configurations that go beyond 
the RPA. In our specific case these will be lplh(8)phonon configurations. 
In the next sections we shall specify these contribution and derive explicit 
expression for A. Here we already point out, that the final Eqs. (2.9) or (2.19) 
will explicitly include both the ph interaction F and (in the propagator) the 
quasiparticle-phonon interaction. In this respect our formulation differs from 
the approaches containing only the quasiparticle-phonon interaction |65| IHB] , 
and we do not use a phonon representation. 

Secondly, the formulas of the extended theory contain also the quanti- 
ties F and that describe the effective interaction and local charge of our 
"refined" quasiparticles. They play the same role as F and Cg in the TFFS 
and will be parameterized in a similar way. In general they should differ nu- 
merically from the quantities F and Cg because the "dangerous" pole terms 
corresponding to the complex configurations, which will be considered ex- 
plicitly in the propagator A, have been removed from F and e^. As discussed 
in section 1.4, however, these differences turn out to be small, as we shall see 
in sections 3 and 4. 

Finally, we have written Eqs. (2.9), (2.20) and (j2.19|) in a new shell model 
representation A, which differs from the representation used in the TFFS. 
These new single-particle wave functions (fx and energies ex that represent 
the new basis {(p\, e\} instead of the phenomenological one {(px, ca}, do not 
contain the phonon mixing, which is included in the mass operator S*^. In a 
sense, the basis corresponds to a Hartree-Fock basis. The mass operator E 
defines the basis functions of the configuration space. The reason for this is 
that in the GF formalism the basis {(fx, ex} is chosen in such a way that it 
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diagonalizes the Greens function (5, see Eq. (2.6): 

Gi2(e)=5i2Gi(e), Gi(e) = -, S ^ +0. (2.21) 

e — ei + ^c^lO 

In this connection the mass operator is nothing else than a generahzed single- 
particle potential. In Eq. (2.21) we introduced the quantity cxi, which is equal 
to +1 for particles and -1 for holes, and is related to the level occupation 
number by cri = 1 — 2ni. The "refined" basis {(fx, ex} can be calculated if 
the old phenomenological basis {(fx, Ca} and the quantity S'^ are known. In 
section 3.1.2 we shall describe our method of calculating {(fx, ex}- 

2.2 Approximation for the mass operator 

and the irreducible ampUtude in the case of 
lplh0phonon configurations 

Let us return to Eq. p.l9|) . which is the starting point of our approach. 
We have shown in section 1.2.3. that by an appropriate choice of the ph 
propagator A and the effective interaction F, that is, in fact, using the ap- 
proximations S = S and U = U, one obtains from this equation Eq. (1.31) 
for the renormalized response functions R, which is the basis of Migdal's 
theory. Therefore our problem of including more complex configurations in 
the generalized renormalized response function defined by Eq. (2.9) requires 
the appropriate choice of the corresponding quantities E*^ and If^, which de- 
scribe the coupling of ph configurations with the more complex ones and 
represent, therefore, the basic input to our extended theory. Depending on 
the approximations made for these quantities, we obtain the various versions 
of extended RPA theories mentioned in the introduction to this chapter: 

1. The Extended Second RPA. If we include in S*^ and If^ complex con- 
figurations in the lowest order of perturbation theory in the effective 
nucleon-nucleon interaction, we obtain the Extended Second RPA of 
Ref. [88]. Here one may or may not consider GSC Such a variant of 
inclusion of "pure" 2p2h configurations within the GF method is briefly 
discussed in ref. jHOj. As it is known that such a perturbation series 
in the nucleon-nucleon interaction converges very slowly, one has to 
consider partial summations that include the interaction to all orders. 
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2. The Phonon Coupling Model. Such a possibility is given by the phonon 
couphng modeL Here the interaction between the ph pairs that build 
up the phonons is included to all orders. Moreover, in the case of closed 
shell nuclei, one can show that the coupling between the particle (or 
hole) with the phonons, the dimensionless parameter g (Eq. (1.44)), is a 
small quantity. For that reason one can restrict oneself to the coupling 
of one phonon to that particle (or hole), which is the g^ approximation 
discussed earlier. In the g^ approximation the expressions for and 
are represented graphically in Fig. 2.1 and have the following analytical 
form: 

^uie) = E T AT (2.22) 

m{a)* m{u) 

^i\34(-, ^, e') = E ^ 7" (2.23) 
e - e' + a[ujm - id) 

Here we have introduced the abbreviation: 

■m(a) _ r m , r m* 

9l2 — 'J(T,+l9l2 + Ca,-lfl'21 ■ 

where g^ is a matrix element of the creation phonon amplitude and 
Um is the phonon energy. Upon substituting Eqs. (|2.22p and (|2.23|) into 
Eqs. (2.1), (2.2) and (2.4), (2.5), (2.7) one obtains the propagator A 
and Eq. (2.9) for the renormalized response function in the ETFFS. 

The simplest propagator A, obtained from Eq. (2.12) in the g^ ap- 
proximation, is given in graphical representation in Fig. 2.2. The cor- 
responding model, which have used this propagator, was suggested in 
refs. pn and realized in refs. IHH IHS] for the Ml resonances 
in magic nuclei ^^O, '^^'^^Ca and ^osp]^ ^^le A representation. A rea- 
sonably good description of the experimental data, including solution 
of the long-standing problem of the isovector Ml resonance in ^ospj-j 
was obtained. In addition, as mentioned in section 1.3.1, the model 
also reproduced the experimentally known 1"*" levels in ^^O and ^'^Ca as 
a result of including the complex lplh(g)phonon configurations in the 
ground state. Such an approach, however, leads to the following prob- 
lem: the propagator A contains in this approximation poles of second 
order in the variable u at the points of the simple poles of the RPA-like 
propagator. The problem is that the second order poles do not have 
the same physical meaning as the simple poles of the exact response 
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function. This can result in distortion of the energy dependence of the 
strength function near these poles. It was shown, however, that for the 
Ml resonances this distortion was negligible I95j. 

The way out of this difficulty is to perform the summation of the g^ 
terms in the propagator. In ref. ^Hj the method of summation was 
developed both for the linearized mass operator (Eq. (2.22)) and for 
a nonlinear case, which uses the full GF G in the mass operator in- 
stead of G. Even so, the numerical difficulties are still rather serious 
in both cases. There exist, however, a simpler method, which treats 
the summation approximately. The main idea in this method is to per- 
form a summation of an infinite series of some of the g^ terms in the 
propagator. 

2.3 Constructing the generalized propagator: the 

method of chronological decoupling of diagrams 

In this section we consider the method of constructing our generalized ph 
propagator Aiu) [77], where the problem of second-order poles of the model 
inS] has been solved by a partial summation of the diagrammatic series 
in the propagator. 

2.3.1 The basic idea. 

Here we present the technical aspects of the method we use to construct the 
generalized ph propagator A{uj). As in section 1.2.2, we denote the set of 
one-particle quantum numbers and the time variable by 1 = {i,ti} and so 
forth. In particular, in the time representation the Green function G given 
by Eq. (2.21) becomes 

G(l, 2) = -i<y^5i2e{aih2)e-'''''\ (2.24) 

where ti2 = ti — ^2 and is the Heaviside (step) function. 

Let us consider an operator in the time representation whose matrix ele- 
ments in the basis {(p\, e\} are given by 

D(12, 34) = 5.„_.,^(ait4i)^(ait32)G(3, 1)G(2, 4). (2.25) 
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Multiplying the right-hand side by the sum 0{aiti2) + 0{ait2i) = 1 and using 
Eq. (2.24) we find 



(2.26) 



D(12,34) = 5.„_.,5i3524e-<^'^*^^+^"^*^^) 

X [0iaiti2)e{aiUi)e{aitsi) + ^(ait2i)^(ait32)^(^it42)] 

Here we have used the easily verified identity (cr = ±1), 

e{atis)e{ah2)e{at23) = e{ati2)0{at23). (2.27) 

We further introduce, in accordance with Eq. ()1.19p and ref . ^H] ? the time 
difference variables ri = ts — ti, r2 = t2 — ^i, = — ^4, so that 

D(12,34) = £)i2,34(ri,r2,r3), (2.28) 

and transform to the energy representation: 

dndT2dTse'^^^^+'^^+'^^^D,2MruT2,T3). (2.29) 

-oo 

Substituting Eq. (2.26) into Eq. (2.29) and using the well-known formulas 



dTe'^^^eiar) = -, a = ±1, 6 ^ +0, (2.30) 

oo Lu + tab 



we find 



/^i2,34(u;,e,e') - 5.„-.,5i3524 + ^ _ + _ + ,^^5) x 

1 1 

e + — ei + iai6 e — 62 + «a2^ 

/^From this, using the definition Eq. (2.21), we obtain the final result 
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-Di2,34(t^, e, e') = _<^2(5i3524cxi(t^ - ei2 + «(Ti5) 
X Gi(e + cu)G2(e)G3(e' + cu)G4(e'). 



(2.32) 



We see from this relation that the dependence of the function 1^12,34(0;, e, e') 
on the variables e and e' is separable. This is essential for the subsequent 
model transformations. We note that the simple product of two Green func- 
tions G(3, 1)G(2, 4) in Eq. (2.13) does not possess this property in the energy 
representation: 

^?2,34(^, e') = -5i352427r5(e - £0^1(6 + uo)G2{e). (2.33) 
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In order to illustrate the effect of this separability, let us consider a func- 
tion -F(12, 34) which, in the time representation, is given by 

F(12,34)= J2 ^^(12, 56)D(56, 78)^^(78, 34), (2.34) 

56,78;t 

where and are different functions whose dependence on the time 
arguments reduces to a dependence on three time difference variables, e.g., 

V^^(12,34) = V,l,,{h,,t2uh^) (2.35) 

and J2i2...:t denotes summation over the one-particle indices of the basis wave 
functions (fi and integration over the time variables. Then, in the energy 
representation, Eq. (2.34) is given by 

^12,56 (^^ e, ei)D56,78(w, Cl, e2)V^^8,3A{^ ^ ^2, e'). 

(2.36) 

Substituting Eq. (2.32) into this, we see that, owing to the separability of 
the energy dependence of D^Qjsi^, ^i, 62), the integrations over the variables 
ei and 62 in Eq. (j2.36p decouple. From the technical point of view this is due 
to the presence of the additional (compared to i?°) factor of two 6 functions 
on the right-hand side of Eq. (2.25) written in the time representation. The 
replacement of by -D in the expressions determining the generalized ph 
propagator is the main element of the method of constructing the propagator 
in our approach. The physical meaning of it will become clear with the 
specification of the quantities entering into equations such as Eqs. (2.34) and 
(2.36) and in terms of Feynmann diagrams. This also provides the name of 
the method of constructing of the generalized propagator, which we refer to 
as the method of chronological decoupling of diagrams (MCDD) [77j . 

2.3.2 Formulation of the method 

Using symbolic notation one can write the solution of Eq. (2.12) in the form 

W = R^ + iR^T^R\ (2.37) 
where the amplitude F*^ satisfies the equation 

= W^ + iW^R^T\ (2.38) 



38 



We introduce the new amplitude defined by the equation 

= W^ + tW%-D)f'', (2.39) 

which is obtained from Eq. (j2.38p by the replacement 

^ -D,W^ -^W". (2.40) 

Here the function D is given by Eq. (2.25) and the amplitude is given by 

W'iU, 34) = W^{12, 34) + I^^°"^P(12, 34), (2.41) 

which differs from Eq. (2.14) by the replacement of — zS'^S'^ by W^"^^. The 
amplitude in Eq. (2.41) plays the same role as the term — zS^S^ in 

Eq. (2.14) and will be defined below. We now replace in Eq. (2.37) the 
amplitude by V^: 

= rO^ ^fiofe^o ^2.42) 
In the energy representation this function defines the desired propagator, 

Ai2m{uj) = j^^ _i?^2,34(^,e)- (2.43) 

Let us discuss the physical meaning of these model transformations. It is 
clear that j^comp contains g^ terms (see also Eq. (2.44), below) and that the 
solution of the integral equation (2.39) (after the substitution of Eqs. (2.22) 
and (2.23)) gives the function R^, and therefore the propagator A, which 
contains an infinite sum of the terms of higher order in the g^ terms. An 
additional — and physical — condition, whose fulfilment must be verified, is 
that this sum must correspond to a subset of the series of Feynmann graphs 
that correspond to the initial quantity R^. In this sense the replacement of 

(Eq. (2.13)) by -D (Eq. (2.25)) is justified by the fact that the quantity 
—D contains 9 functions, which are projection operators in the space of time 
variables; that is, the graphs corresponding to different combinations of argu- 
ments of the function — D(12,34) belong to the set of graphs corresponding 
to the function R^{\2^ 34). Physically, however, it is very important that the 
replacement of RP hy —D eliminates from the diagrammatic expansion of i?^ 
the terms corresponding to processes in which configurations more complex 
than lplh(8>phonon are excited while leaving most of the lplh®phonon con- 
figurations. In our model with IplhCSphonon configurations the discarded 
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graphs include those in which a time cut through two or more phonon hues 
is possible in the time representation. As an example, in Figs. 2.3a, 2.3b and 
2.3c we show three graphs of order g^ that are excluded from the expansion, 
and in Figs. 2.3d, 2.3e and 2.3f, we show three similar graphs that remain in 
the diagrammatic expansion of R^. 

The sum of the contributions of higher order in g^, which is contained in 
the function R^, corresponds to the sum of chains of graphs similar to those 
shown in Figs. 2.3d, 2.3e and 2.3f. The contributions of the graphs shown in 
Figs. 2.3a, 2.3b and 2.3c are excluded due simply to the 9 functions in the 
definition of Eq. (2.25). 

The statements that in the method under consideration all the contribu- 
tions of configurations more complex than lplh®phonon are excluded, while 
all lplh®phonon contributions are included, are completely valid only when 
GSC are ignored. The role played by their effect will be discussed below. 
These effects correspond to the so-called backward-going diagrams which, 
on the one hand, take us beyond the lplh®phonon (or 2p2h) approxima- 
tion (see Fig. 2.4) but, on the other hand, not all the lplh®phonon con- 
tributions related to them (those of order g^ and above) are included in the 
MCDD. Since g^ is a small parameter, however, these neglected contributions 
of higher order are small. Our task in constructing the generalized propaga- 
tor is to avoid the second order poles, that is, to perform approximately an 
infinite summation of the graphs of the type shown in Figs. 2.3d, 2.3e and 
2.3f and, simultaneously, to take into account the remaining graphs — at least 
within the accuracy of g^. 

The last formula, which is necessary to obtain our generalized propagator, 
is an expression for W'^""^^ in Eq. (2.41). The role of this quantity is to remove 
the multiple counting of graphs with self-energy insertions S*^. This problem 
therefore arises only for the terms of order g^ and above and for the inclusion 
of the GSC induced by lplh®phonon configurations because, if there are no 
backward-going graphs like those shown in Fig. 2.4, the multiple counting is 
eliminated by the 6 functions in the definition of the function D (Eq. (2.25)), 
and in that case W™™^ = 0. The expression for W'^°^^ was found from the 
condition of cancellation of the double-counting of self-energy contributions 
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obtained after the first iteration of Eq. (2.39) [77] : 

El'2'3'4';i E5'6'7'8';t ^"1(1', 1)G(5', l')G-l(2, 2')G{2', &) X 

S^(7', 5')S^(6', 8')G{3', 7')G'\3, 3')G{8', 4')G-\4', 4)e(l'2', 3'4'), 

(2.44) 

where 9(1'2',3'4') = 6'((Ji/t4/i/)^(cri/t5/8/)6'(cri/t3'2')^(c'i'i6'7')- This expression 
was used in our calculations (the explicit form of it is given in the Appendix). 
An improved form of this term was discussed in ref. [20] • AH these effects, 
however, are of order of g^, so that the difference is small in our case. 

2.3.3 The generalized propagator in the energy representation 
and its properties 

Equations (2.25), (2.39), (2.41), (2.42), (2.43) and (2.44) completely deter- 
mine the generalized ph propagator A{uj) in the MCCD. Since most of them 
are written in the time representation it is necessary to transform to the 
energy representation. After some algebraic effort we obtain 



^12,34(u;) = 

E56,78[^15526 + QSts*] (^)]46J^ (^) [573^84 + Q^'s^M] + PittK^) 



(2.45) 



where 



Qitui^) = -5<.„.2<^.3,-<.4^12,34(^)^3(^ - 634) (2.46) 

Q^uili^) = -5.„-.,5.3,.4^rj,34(^)^i(^ - eu) (2.47) 
PlttlH = 5.„..5.3,<x4^S,34(^) (2.48) 



^Sl34(^) = - -^G^ie + uj)G2ie)U!,,,,iu, e, e')G,ie' + 0^)^4(6')- 

J^^^''^'^' + ^)G2{e)^l,{e)G,{e) + 
^24^3(6 + uj)El,{e + uj)G,{e + uj)G2{e)], 

(2.49) 
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with ei2 = ei — e2 and the GF Gi{e) given by Eq. (2.21). 

The only equation that must be solved to find the propagator A{uj) is the 
equation for the ph-ph part ^sejs of the propagator: 

^12,34(^) = ^12,34(W) - E ^2,56(^)'^'56,78(^)48,34(w), (2-50) 

56,78 

where the propagator 

^12,34(t^) = -5l3'^240-l5ai -aa/ " ^12) (2.51) 

is the "refined" RPA propagator and 

$12,34(^) = <i'l2,34(^) + (2.52) 
^12,34(C^) = -5^1 -a2<5a:i -a4Cri(^^ " £12)^12,34 (^)o'3(t^ " £34)- (2.53) 

The expression for ^'^°"^p[uj) is given in the Appendix. 

Properties of the generalized propagator A{uj) have already been dis- 
cussed in detail in ref. j^ni- Here we enumerate only the most important 
ones. Let us first introduce for convenience a new terminology. Let us rep- 
resent our propagator in the form 

A{uj) = )(cu) + A^~+\uj) + A(+")(cu) + A^++\uj) 

A(-+\uj) = A(-~\uj)Q^~+\uj), A(+-)(cu) = g(+-)(cu)A(— )(cu), (2.54) 

A^++){uj) = g(+-)(cu)A(--)g(-+)(cu) + P(++)(cu) 

and refer to the term A^ ^ as the unlike component (e.g., the ph-ph com- 
ponent), and to the remaining three terms as the assocmted components. 

First we consider the question of the second order poles, which is the main 
technical problem that the MCDD is designed to solve. It is obvious, in fact, 
that ^(a;) does not contain second-order poles at the points u = ±(ep — ih), 
which are simple poles of the RPA propagator A{u!). Moreover it follows from 
Eq. (2.50) that the function A^ \ and consequently A{u;) too, is regular at 
these points, as would be expected if the matrix [^"^(0;) + $(ci;)] is invertible 
at these values of u, as is true except for accidental cases. One can also prove 
PU] that: (i) in solving the problem of constructing the MCCD propagator 
in the complete one-particle basis, all its components are regular at singular 
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points of the function A^^^cu) (Eq. (2.49)), and (ii) when the full one-particle 
basis {01, ei} is used, all the poles of the propagator A{u) coincide with the 
poles of its unlike component K^)- 

As for the poles of the function ^(cj) itself, the number of poles of 
(llj) can be greater than the total number of poles of the functions A{uj) 
and ^{u) due to splitting of the poles of $(a;) in the complete or partial lifting 
of their degeneracy in the solution of Eq. (2.50). These effects are of order g^. 

Let us also point out a role of the term Q^~^~'^A^ ^Q*^ in Eq.(2.54) 
for the MCDD propagator, although formally this term is of order of g^. 
The analysis of its poles shows that they coincide with the poles of the 
unlike component A^ ^ and can ensure the removal of the singularities of 
the function A^~^^\ which coincide with the poles of p(++) of Eq. (2.54). 
This absence of "extra" poles provides an indirect confirmation of our model 
approximation. 

Using Eqs. (2.45)-(2.49) or formulas from the Appendix, one obtains for 
the MCDD propagator 

E ^11,34(^0') = E ^12,33(^) = 0. (2.55) 
1 3 

It follows from Eq. (2.9) that analogous identities are valid for the renormal- 
ized response function R{uj). We then find that the relation 

Epii = (2-56) 
1 

holds in our approach. This means that an external field does not change the 
total number of particles; that is, we have the particle number conservation 
law. Were the problem to be solved exactly, this would be obvious from very 
general conditions. As in many other models, however, the validity of the 
law is not obvious a priori. This implies that the graphs with self-energy 
insertions E*^ and phonons in the cross ph channel If^ must be taken into 
account simultaneously. 

2.3.4 Ground state correlations induced by lplh(g)phonon config- 
urations 

As discussed earlier, in the GF graphical language any GSC correspond to 
the so-called backward-going diagrams. In Fig. 2.5 we show some examples of 
the graphs of order g^ that correspond to the GSC caused by lplh®phonon 
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configurations (GSCphon)- In the approach under discussion there are two 
types of GSCphon effects: those wliose inclusion, as in the RPA GSC case, 
affect only the location and the magnitude of the residues of our ph propaga- 
tor A(ct;) and those of new type that can give additional poles of the function 
^4(0;). In the latter case our approach is qualitatively different from the RPA. 
Some of the corresponding graphs are shown in Figs. 2.5e-h. These new poles 
are caused by the associated components of the MCDD propagator, which 
are directly connected with GSCphon effects. The same components can also 
give GSCphon effects that are not connected with the new poles and, in this 
sense, are analogous to the RPA GSC. 

Another important feature of the associated components of the MCDD 
propagator is that the solution of Eqs. (|2.9|) or (2.19) using this propaga- 
tor also contains associated components. Therefore — and contrary to the 
RPA case, according to the spectral expansion of the response function, 
non-zero components of the transition densities Ppp, and p^^/ appear, where 
P12 =< ?T.|a^a2|0 >. These effects are caused by the GSCphon corresponding 
to the above-mentioned new poles. Physically, the corresponding graphs de- 
scribe processes of creation from the vacuum or annihilation of two ph pairs 
simultaneously, i.e, 2p2h configurations. In other words, the external field 
can directly excite 2p2h configuration from the correlated ground state in 
such a way that a particle of the 2p2h configuration makes a transition from 
one state to another {pp' transition) or a hole of this configuration makes 
a transition to another hole state {hh' transition). This means that new 
components, Ppp, or p^^,, appear in the spectral expansion of the response 
function in addition to the usual and p^^ components. The inclusion 
of these effects leads to a change of the sum rules for the moments of the 
strength function. Within the ESRPA, that is, for the case of "pure" 2p2h 
configurations, this question has been treated analytically in ref. j^ni and in 
ref. jSni; where calculations for the Ml and Gamov- Teller resonances in '^^Ca 
have been reported. Our calculations (see section 3), where perturbation 
theory on the ph interaction was not used, show the quantitative importance 
of the inclusion of the associated components of the propagator. 

We can already see from this brief discussion of the GSCphon effects that 
they manifest themselves in a considerably more complicated way than the 
RPA GSC effects. This question therefore merits further study. 
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2.4 Inclusion of the single-particle continuum 

In the coordinate representation the equation for the density matrix is ob- 
tained from Eq. (2.19): 

p(r, uj) = -J A{r, r', uj)e,V^{r', uj)d^r'- 

/ A(r, ri, w)F(ri, r2)p(r2, a;)dViC?V2. 

The present equation for the change of the density matrix is formally 
identical with the corresponding equation (1.35) for the case of the standard 
TFFS. However, the quantities which enter {cg, F} and the propagator A 
of the extended theory are defined differently compared with the quantities 
{eg, F} and A of the standard theory. The propagator of the ETFFS differs 
the most from the standard RPA propagator. It consists of two parts: the 
RPA-like part and the part that contains much more complicated physics 
caused by lplh®phonon configurations, including GSCphon- 

As discussed for the case of the standard TFFS (section 1.2.4), use of the 
coordinate representation allows accounting for the single-particle continuum 
exactly. The effects associated with resonance decay into the continuum (the 
width r I), that is, with nucleon emission from the nucleus, are thereby 
automatically included in the width. 

The systematic use of this technique for the case of complex configura- 
tions leads to great computational difficulties. (See, for example, the general 
expression for A{r,r',Lj) in the g^ approximation obtained in ref. [HI]-) At 
this stage we therefore included the continuum only in the RPA part of the 
propagator, using the idea of the so-called combined (r. A) representation 
that was developed within the TFFS for nuclei with pairing [^Tj. In our 
ETFFS the expression for the propagator A is given by 

A{r,r',uj)=AZt{r,r\uj) + 

Ei234(A234(^) - A^it{um{r)Ur)Ur')^lW) , ^ ' 

where A^f is the refined RPA propagator in which the single-particle con- 
tinuum is taken into account exactly. In the following applications the sum- 
mation in Eq. (|2.58p is performed over two shells above and all shells below 
the Fermi level. If one would sum over all configurations, the terms with 
j^RPA ()2.58|) would cancel each other and we would have the exact 

expression for A{r, r', u). It is necessary to subtract the Af^^f term in order 
to avoid double counting as the full propagator A1234 contains the RPA part 
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already. Thus, in this method the single-particle continuum is taken into 
account completely only in the RPA part of the generalized propagator. Our 
calculations show that this approximation is satisfactory. The expression for 
the propagator A1234 in the "refined" single-particle basis {(p\,e\}, is given 
explicitly in the Appendix. 

The distribution of the transition strength is then given by 

S{u, A) = -- Im / d^r[egV%r, u)]*p{r, u + lA) , (2.59) 

from which one can easily obtain the transition probabilities and the EWSR 
summed in an energy interval. 

Equations (2.57), (2.58) and (2.59) are the basic relations used in the 
calculations reported below. 



3 Application to giant resonances 
3.1 Scheme of the calculations 

In the previous section the basic equation (2.57) of the ETFFS for the case 
of closed shell nuclei has been formulated in coordinate space for the change 
of the density matrix in an external field V^{r,uj) with the energy u. The 
strength function (Eq. (2.59)) gives the energy distribution of the excitation 
strength under consideration: 

g(cu. A) = "^^[^^^ = -- Im n(cu + zA), (3.1) 

where 

U{uj + tA) = J d^r[eqV\r)]*p{r, lu + lA) (3.2) 

is the polarizability propagator and A is a smearing parameter. By using this 
we take into account phenomenologically those complex configurations that 
are not explicitly treated and we simulate the finite experimental resolution. 
The use of a sufficiently large A also greatly reduces numerical difficulties 
in the calculations. From the solutions of this equation one obtains the 
transition probabilities and the EWSR, summed over the energy interval 
[El, E2], —e.g., the hnear EWSR: 

Sl= EMEL) T = T^fdoo ooUioo). (3.3) 
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Here the integration contour in the complex plain intersects the real energy 
axis at the points Ei and E2. 



3.1.1 Electric sum rules 

In the following we discuss electromagnetically induced EO, El and E2 ex- 
citations. In order to avoid possible confusion we will define all relevant 
quantities exphcitly. The electric operators are given as: 

z 

QLM = e^rfyiM(a), L>2 (3.4) 

i=l 

Z 

Qm = eY.r} (3.5) 

i=l 

eN ^ eZ ^ 

Qui = -j-Y.'^i'^^Mi^i) - -j-Y.'^i'^^Mi^i), (3.6) 
1=1 1=1 

where the electric isovector dipole operator contains the well-known kine- 
matic charges Cp = Ne/A, e„ = —Ze/A. The linear (full) energy-weighted 
sum rules for these operators 

Sl = EWSR = Y.{Ek - Eo)Bk{EL) j 

k 

are then given by |117j 

Sl = ^^L{2L + ifZ < r^^-^ >p, L>2, (3.7) 



So = Z < >p, (3.8) 



nip 



Snip A 

where the radial average is usually taken over the proton distribution in the 
ground state. The electric dipole EWSR, Eq. (3.9), corresponds to the well- 
known Thomas-Reiche-Kuhn sum rules for photoabsorption cross section. 



= J a{E)dE 



2iT^e^hNZ 
mc A 
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The sums 5*2 and 5*0 include both the isoscalar (AT = 0) and isovector 
(AT = 1) contributions, where T is the isospin of the nucleus. These contri- 
butions may be separated for the EO and EL (L> 2) transitions by means of 
the additional physical approximation < > = (A^ — Z^j A |117j so that 



Qlm — Q\m + Q\mi 



where for L > 2 



and 



QLM = ^ErflLM(l^.) (3.10) 

1=1 



and for L = 



eN ^ eZ ^ 

QIm = -^T.r^'YLMm ~^Y.^lY,Mm (3.11) 
^ i=l ^ i=l 

Qlo = "-^j:rl (3.12) 

eN ^ eZ ^ 
Qlo = -^Er^--fErl (3.13) 

i=l i=l 

The effective charges in Eqs. (3.11) and (3.13) coincide with the kinematic 
charges for isovector El excitations, but they have another physical origin. 
These operators give the following isoscalar and isovector electric EWSR 

STTirip A 

for L > 2 and 

^<r^>. (3.16) 
, 2h^e^NZ . 

for L = 0. 

We have calculated our strength distributions S{E) with the electric op- 
erators (3.4), (3.5) and (3.6) for the full S{E), (3.10) and (3.12) for the 
isoscalar S{E) and (3.11) and (3.13) for the isovector S{E). 
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3.1.2 Input quantities for the ETFFS 

As discussed earlier, as in tlie standard TFFS, the initial numerical input to 
the ETFFS are two sets of phenomenological parameters that describe (i) 
the Woods-Saxon single-particle potential and (ii) the effective interaction of 
the Landau-Migdal type. 

In all of the approaches that include complex configurations in this way, 
two problems concerning double counting of the complex configurations (in 
our case, lplh®phonon) arise. They concern the necessity to extract the 
contribution of the complex configurations that are treated explicitly from 
the "old" phenomenological quantities, that is,(i) from the mean field and 
(ii) from the effective interaction and local charges of quasiparticles. 

As far as the procedure for the mean field is concerned, this "refining" 
from our phonons goes as follows: To obtain the new ("refined") single- 
particle basis {ex, ipx} we must subtract the contribution of the phonon mix- 
ing from the "old" phenomenological Woods-Saxon energies PTI I92j . 



where S*^ is the energy-dependent part of the mass operator in Eq. (2.1). In 
our g^ approximation it is defined in Eq. (2.22). More detailed expressions 
are given in the Appendix. 

As for the effective interaction and local charges of quasiparticles, in all 
our equations the new effective interaction F and local charges Cg enter and 
play the same role as the corresponding quantities F and of the TFFS. 
By definition they should not contain the complex configurations that are 
considered explicitly and, therefore, the parameters describing in this ap- 
proximation F and have to be determined, in principle, from experiment. 
However, as discussed in section 1.4, we use only a small number of low-lying 
phonons, which give the main contribution to the giant resonances charac- 
teristics under consideration. The corrections to the quantities F and Cg due 
to these phonons must be irregular and of a long-range character. Therefore 
one can hope that these non-local corrections of the local quantities F and 
Cg are not very important and that we can use the old parameters of F and 
Cg instead of the new ones. A detailed investigation of the Landau-Migdal 
parameters shows indeed that such changes, except in the case of fex, are 
negligibly small. 

Thus, although additional configurations, are treated in the ETFFS ex- 
plicitly, it is possible to use the the parameters of the TFFS approach. This 
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means that all but one, of the input quantities of the ETFFS are, at least to 
a good approximation, already known. 

3.1.3 Numerical details 

For the reasons mentioned above, we use in the following application of the 
ETFFS the same Landau-Migdal interaction as we had used previously in 
the standart TFFS, where we also investigated low- and high-lying collective 
states. The parameters are given in Eq. (1.42). The only exception is the 
parameter f^x- Because of the non-self-consistency of the Landau-Migdal 
approach, fex should be chosen in such a way that the spurious 1~ state lays 
at zero energy. In our calculations we used for the density interpolation in 
Eq. (1.11) the theoretical ground state density po{r) given by Eq. (1.43). In 
our opinion, such a choice of po{r) instead of the usual Woods-Saxon form — 
which we have used starting with ref. jHS] — makes our calculations more 
consistent. (See also the discussion of previous values of fex in ref. pUj.) 

For the medium mass nuclei ^°Ca, ^^Ca and ^^Ni, using Eq. (1.43) gives 
noticeably different results (mainly for the IS EO resonances) from those ob- 
tained using the Woods-Saxon form of po{r) as compared with the difference 
for ^°^Pb. In particular, this change in the numerical scheme has shown 
that the effects of GSCphon there were numerically overestimated by us |137j 
although their qualitative features were left unchanged |H3] . 

The phonon characteristics, that is, their energies and transition proba- 
bilities, which have been calculated within the standard TFFS j2] with the 
same parameters (Eq. (1.42)), are given in Table 3.1 and in ref. |137j . The 
criteria for the selection of these phonons were that they should have the 
strongest coupling parameters — and therefore the largest transition proba- 
bilities, and that the energies of the corresponding configurations should be 
in the energy region of the giant resonance. 

The influence of the different number of phonons under consideration 
was investigated for ^°^Pb. In ref. j^H] it was shown that the role played 
by the "unnatural parity" states (e.g., 2j",3]'',4j~ and 5f) is very small. A 
similar result for a much larger number of unnatural parity phonons has been 
obtained within a simpler model [HHl- In refs. (HHIIZHI the contribution of the 
four 2"*" phonons corresponding to the isoscalar E2 resonance at an energy 
of about 10 MeV in ^°^Pb has also been checked. It turned out that their 
contribution was only noticeable in the high-energy tail of the El resonance. 

The single-particle levels and wave functions were calculated with the 
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standard Woods-Saxon potential |l(JUj . In order to get good agreement with 
the experimental single-particle energies, the well depth of the central part of 
the potential was adjusted by changing the depth parameter by less than 5 
percent. The single-particle energies thus obtained are given in refs. |8UllH3] . 

Numerical values of all the EWSR used in the calculations are presented in 
Table 3.2. The quantities < >p were obtained using the proton density dis- 
tribution in the ground state, Eq. (1.43), calculated from the single-particle 
scheme described above. 

The mean energies and dispersion D are determined by 

mo mi mo 

(3.19) 

where the energy moments rrik corresponding to the energy interval AE = 
Emax — Emin havc been calculated using 

mk= 1^^^^ dE E^S{E). (3.20) 

3.2 Electric and magnetic resonances in magic nuclei 
3.2.1 Photoabsorption cross sections; El excitations. 

The El absorption cross section has been calculated using the formula 

(Tsi(cu) = 4.022 cj5si(cj), (3.21) 

where uj is taken in MeV, 5* in fm^-MeV~^ and a in mb. 

The results of these calculations in ^°Ca, ^^Ca and ^osp]-, ^^^^ ^j^g com- 
parisons with experiment are presented in Figs. 3.1-3.3 and Table 3.3 [75] . 
The label lplh-|-continuum in Table 3.3 denotes the calculations done within 
the CRPA (or CTFFS) and the label Iplh-h2p2h denotes those without the 
continuum but with the smearing parameter A 7^ 0. In order to obtain the 
integral characteristics of the resonances we used a Lorentz function to ap- 
proximate the resonance curves, as is usually done in analyzing experiments. 
The parameters of this function — the mean energy E, the maximum value 
of the cross section amax and the resonance width F — were found from the 
condition that the three lowest energy moments in Eq. (3.19) {k = 0,1,2) 
coincide for the exact and approximate resonance curves. The characteristics 
of the fitted lorentzians — the mean energy E, the width F and the integral 
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cross sections ctq — were obtained for the same experimental energy interval 
and are compared with the corresponding experimental values in Table 3.3. 
For details see ref. [75] . 

In the case of CRPA we see that the mean energy values E of the El 
resonances and the widths do not agree with experiment. The energies are 
too low by 1 to 2 MeV, where as the widths are off by more than a factor of 
two. The inclusion of the lplh(8)phonon configurations noticeably improves 
the agreement with experiment as compared with the CRPA. As can be seen 
from Figs. 3.1-3.3, the complex correlations lead to a large reduction of the 
maxima of the photoabsorption cross sections and to a better description of 
the El resonance tails. 

The final results (denoted in Table 3.3 as lplh+2p2h-|-continuum) show 
that inclusion of the continuum changes the values of E by 0.6-0.8 MeV 
in the required direction as compared with the lplh-|-2p2h approximation, 
except for the case of ^°Ca. The theoretical mean energies for ^^Ca and ^osp]-, 
are now in good agreement with the data. 

The most important result of these calculations, however, is the now 
satisfactory explanation of El widths. To be precise, they are explained to an 
accuracy of the value of the smearing parameter used, which is smaller than 
the widths. The two new (as compared with the standard RPA) elements, 
namely the lplh(8)phonon configurations and the single-particle continuum, 
increase the values of P by about a factor of two. As should be expected, 
the contribution of the continuum to P is small for the heavy nucleus ^ospb 
and somewhat larger for the lighter ones. Specifically, its contribution to the 
calculated total width is 14 percent for ^°Ca and 28 percent for ^®Ca, but only 
7 percent for ^''^Pb. Now we underestimate Pexp by only 14 percent for ^''Ca, 
11 percent for ^^Ca and 4 percent for ^°^Pb. These small differences may be 
connected with the approximate treatment of the single-particle continuum. 

In Figs. 3.1-3.3 the theoretical curves were approximated by lorentzians 
in the energy interval in which there are experimental points and, corre- 
spondingly, the integral characteristics of Table 3.3 were calculated for the 
same interval. In order to study the EWSR questions and the role of the 
GSCphon it is necessary, at least at first, to increase the energy interval un- 
der consideration. Therefore, in our later calculations [22] for ^°Ca, ^^Ca and 
^^Ni, considerably larger energy intervals were studied and the corresponding 
theoretical curves are presented there. Comparing these results with experi- 
mental curves for ^°Ca [103j and ^^Ca |104j . one can conclude that we have 
also obtained agreement with experiment for mean energies, total widths and 
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maximum values of photoabsorption cross sections, although the agreement 
is better for the previous analysis, as should be expected. This agreement is 
mainly obtained by inclusion of lplh®phonon configurations, but the role of 
our GSCphon is noticeable, for example, in the integral quantity El EWSR 
(see Table 3.4. and section 3.2.4.). We have also calculated the depletion 
of EWSR (Jo in the experimentally studied intervals 10.0-32.0 MeV (ctq^^ = 
637.7 mb-MeV) for ^°Ca and 11.0-27.5 MeV (ctq''^ = 836.6 mb-MeV) for ^^Ca 
and obtained reasonable agreement with experiment: 599.2 mb-MeV (99.9 
percent of a^) and 635.3 mb-MeV (90.8 percent of a^), respectively. 

The integral characteristics calculated for the large intervals are given 
in Table 3.4. The depletion of the corresponding EWSR (in percentages of 
Si, Eq. (3.9)) in Table 3.4 were obtained using the values of the reference 
EWSR from Table 3.2. We have obtained 97-101 percent of the El EWSR for 
the RPA well as for the RPA+lplh®phonon configurations (without 

GSCphon)- The case with GSCphon will be discussed in section 3.2.4. 

3.2.2 EO and E2 resonances 

The isoscalar electric monopole (IS EO) giant resonance in nuclei is a unique 
source of information on the compressibility and equation of state of nuclei. 
The extrapolation to nuclear matter and neutron stars requires that the 
energy of the resonance be known over a wide range of the mass number A and 
for very different numbers of protons and neutrons. Here we discuss ETFFS 
results simultaneously for the EO and E2 resonances because these resonances 
have rather similar excitation energies and overlap to a large extent. The 
results of isoscalar (IS) and also isovector (IV) EO and E2 resonances are 
given in Figs. 3.4, 3.5, Tables 3.5, 3.6 for ^°Ca, ^^Ca and ^^Ni and in Figs. 3.6, 
3.7 and Table 3.7 for ^ospb. 

EO and E2 resonances in ^°Ca and ^^Ca 

The main part of the IS EO-resonance in ^°Ca — 65 percent of Sq as taken 
from Table 3.2, is in the 11.0-23.0 MeV interval (see Fig. 3.5). In the 6.0- 
30.0 MeV interval we have 85.5 percent.. The comparison with experiment 
for isoscalar EO and E2 excitations^ in ^°Ca has been made in ref. [H2]- R 

"'^In our comparisons with experiments there is an uncertainty because in the experimen- 
tal works, as a rule, there is no exact information about numerical values of the reference 
EWSR used. Therefore there may be a discrepancy due to the different values of < >p 
used here and in the experimental articles. 
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was shown there that in order to explain the electron scattering experiments 
where both of the resonances (EO and E2) are excited, it was necessary 
to take into account GSCphon for both of them. This can be seen in Fig. 3.4 
for the full electromagnetic EO and E2 strengths in ^°Ca and also in Fig. 3.5. 
See also the discussion in section 4. 

For the mean energies of the isoscalar EO resonance in ^°Ca we obtained 
Ef]^ = 18.4 MeV (Table 3.5), compared with the experimental value of -Ei 
= 18.89 ± 0.11 MeV |128j . Our theoretical results agree with experiment 
reasonably well not only for the integral characteristics, but we also for gross 
structure of the isoscalar EO strength function in Fig. 3.5 and for of the cross 
section in refs. |128t 1129] . as we shall see in section 4. 

It was first observed in ref. |109j that the full isoscalar E2 strength in 
^°Ca is located in the 0-22 Mev region and is divided into approximately 
equal parts around 13.5 Mev and 18 Mev. This splitting was confirmed 
later mUllIII]. Our calculations give 73.8 percent of for the 0-22 MeV 
region, 18.3 percent for the 12-15 MeV region and 27.9 percent in the broader 
10.5-16.5 MeV interval. For the 16.5-19.5 MeV and 15-21 MeV regions, we 
obtain 25 percent and 44.6 percent, respectively. It is therefore difficult to 
speak about the egiia/ distribution between the regions around 13.5 Mev and 
18 MeV, but the trend is the same and corresponds roughly to experiment. 
The splitting can be seen in Fig. 3.5. In agreement with the experiment in 
ref. |lllj . we obtained also the third maximum at 12 MeV. 

As can be seen from Fig. 3.5, the isovector EO as well as isovector E2 reso- 
nances are spread out over larger regions than the corresponding isoscalar EO 
and E2 resonances, for one obvious reason: the attractive isoscalar interac- 
tion shifts the isoscalar strength down and that reduces the (escape) width; 
the repulsive isovector interaction has the opposite effect. There are also 
noticeable low-lying tails of isovector strength in the regions of the isoscalar 
EO and E2 resonances. (Table 6 of ref. [HH] contains numerical results for 
all the nuclei under consideration.) The total isovector E0+E2 contributions 
to the main region of the isoscalar EO resonance are about one-half for ^'^Ca 
and ^^Ni and about one-fourth for ^^Ca. These results may be important for 
electron scattering on these nuclei; at least a similar effect has been obtained 
in ^^Si(e,e'a) coincidence scattering |118j . 

Recently, the comparison of our calculations for the low-energy isoscalar 
E2 strength in ^°Ca (in the 10-18 MeV interval) and ^^Ca (in the 11-14 MeV 
interval) with the strength that was extracted from {p,p'x) reactions (x = ao, 
Pq and TT-o) has been performed in ref. |138j . In order to do this the authors 
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multiplied the calculated E2 strength by the ground state branching ratios 
deduced from statistical model calculations. They obtained good overall 
agreement with our calculations, which probably also supports the need to 
take the GSCphon into account. On the other hand, a very large percentage of 
the E2 EWSR observed in both ^°Ca and ^^Ca [Hn] between 6 and 12 MeV 
(about 40 percent and 23 percent, respectively), which means in fact a large 
shift of these E2 resonances — especially in ^°Ca — to the low-energy region, 
is in strong conflict with the results of our calculations shown in Fig. 3.5. In 
these experiments a heavy ion *^Kr beam of 60 MeV/nucleon was used to 
excite the Ca isotopes. 

EO and E2 resonances in ^ospb 

The EO and E2 hadronic response functions in ^°^Pb have been calculated 
here for both the IS and IV resonances within the ETFFS and for large energy 
intervals (see Figs. 3.6 and 3.7 and Table 3.7). The parameter f^x = ^2.2 
was obtained by fitting to the energy of the 2^ level. The difference from the 
value of -1.9 obtained in the CTFFS calculations (section 1.2.5) is caused 
by the complex configurations that have been included here. This small 
difference of the parameters shows the numerical effect of inclusion of our 
complex configurations for this low-lying level in ^°^Pb. 

It is natural (and in accord with the experiment) that, as one can see 
from Table 3.7, the mean energy values and depletion of the EWSR depend 
rather sensitively on the energy interval used in the calculation. In most cases 
we obtained reasonable agreement with the experimental data available for 
these integral characteristics if the energy intervals under consideration are 
comparable. (A recent summary of the experimental data on the EO and E2 
resonances in ^osp]^ jg presented in ref. jHO]-) As for the cases of ^°Ca and 
^^Ca, there is a noticeable amount of the IS strength in the main region of 
the IV strength and vice-versa for both EO and E2 resonances. 

Comparison of the E2 results shown in Fig. 3.7 and Table 3.7 with those of 
the CTFFS in Figs. 1.3, 1.4 and Table 1.1, calculated with the same smearing 
parameter, A = 250 keV, shows the role of lplh(8)phonon configurations: (i) 
its inclusion gives the resonance widths that are numerically similar to the 
observed values of 3.1 ± 0.3 MeV [33] and 5 ± 0.5 MeV [27] for the IS E2 
and IV E2, respectively and (ii) as for the isovector El, the mean energy of 
the IV E2 resonance is changed in the desired direction (see Table 3.7). 

In the calculations with the smearing parameter A = 250 keV shown 
in Fig. 3.7, a gross structure of the IS E2 hadronic response was obtained 
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that is absent, of course, in the calculations with A = 800 keV, shown also 
in Fig. 3.7. Thus one can see a natural and noticeable difference between 
these two calculations that could be checked in hadron experiments. The 
calculated fine structure of the electromagnetic response function is given in 
Fig. 3.9 and discussed below. 

For the EO IS and IV resonances in ^osp]^ calculated with A = 250 keV, 
shown in Fig. 3.6 and Table 3.6, we obtained approximately the same results 
as for the E2 greement of the integral characteristics with experiment, 

as well as a gross structure that probably has not been observed. 

On the fine structure of El and E2 resonances 

Generally speaking, the description of the giant resonance fine structure ^ 
is a natural step once the single-particle continuum and some complex con- 
figurations have been taken into account. A principal motivation for such a 
step is the rapidly improving experimental resolution. The phenomenologi- 
cal smearing parameter used in our and many other calculations in practice 
unifies (or simulates, to be more precise) two quite different effects: realis- 
tic experimental resolution and complex configurations that are not treated 
explicitly in the approach under consideration. 

To exclude the influence of smearing and to imitate, in a sense, the results 
of future El experiments with very good resolution, we have repeated the 
calculations for the El resonance in ^osp]^ with A = 10 keV. The results 
are shown in Fig. 3.8. Of course, it is unrealistic to expect to observe the 
fine structure obtained on the high-lying slope of the El resonance, but the 
one on the low-lying slope could probably be observed if the "strongest" 
complex configurations were chosen correctly. A more detailed discussion of 
the low-energy strength is presented in section 3.2.5. 

The results of the calculations of the E2 resonance fine structure are pre- 
sented in Fig. 3.9 for the electromagnetic IS E2 strength function in ^°^Pb. 
They improve our earlier calculations jHSl in the following way: (i) as in the 
previous calculations under discussion, our G'^Gphon have been taken into 
account; (ii) as in Fig. 3.7, the value of the fex = ^2.2 has been adjusted to 
obtain the energy of the 2f level and (iii) as mentioned above, Eq. (1.43) for 
the nuclear density in the ground state has been used. The smearing param- 
eter was taken as A = 40 keV so as to be comparable with the experimental 

^For definiteness we will refer to an observed structure as fine if the experimental 
resolution is less than 100 keV and gross if the resolution is more than 100 keV. 
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resolution in the (e, e') and {p,p') experiments, results of which are presented 
in Fig. 3.9. 

As can be seen from Fig. 3.9 and corresponding Fig. 1 in ref. (HHl; we ob- 
tained here better agreement with experiment than in ref. [Sn] • In that work it 
was shown that a large part of the fine structure under consideration was due 
to complex configurations. The same result can be clearly seen in Fig. 3.9. 
All of this confirms the decisive role of our complex lplh(8)phonon config- 
urations in the formation of the fine structure — at least for this resonance. 
On the whole, however, the situation is very involved, due in particular to 
the fact that the escape width depends strongly on the excitation energy 
and that there are very many non-collective configurations. If we do not yet 
include additional complex configurations (which is, of course, quite possible 
to do within the GF formalism using the general prescription described in 
section 2.1), then one of the first steps towards clarifying the situation is prob- 
ably to take into account the energy dependence of the smearing parameter 
and/or to use a "refined" optical potential that effectively takes into account 
the complex configurations that we do not consider explicitly [HZl El] • 

In any case, the explanation of the giant resonance fine structure is a 
challenge for any microscopic theory. It is clear that after improving the 
experimental resolution it will be possible to observe more complex configu- 
rations in the fine structure. In this case it will be necessary to improve the 
ETFFS by the addition of new configurations, which should be done follow- 
ing the methods of construction of the ETFFS described in section 2. One 
can also expect successes in this direction for the low-lying structures in the 
nucleon separation energy region (see section 3.2.5). 

3.2.3 Ml resonances 

There is a rich history of study of Ml excitations in nuclei |12Uj . The calcu- 
lations within the standard TFFS with the Landau-Migdal interaction were 
able to reproduce the excitation energies of the isoscalar and isovector Ml 
excitations with the universal values of spin interaction parameters g and g' 
defined in Eq. (1.42) [Ml 1121112111021 • 

The observed total transition strengths were also reproduced in these 
calculations because the standard TFFS contains the universal spin local 
charges determined from [2]: 
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where 7^ = 2.79/io, 7" = — 1.91/io, /io = ^ and 

e = c = 0.1, ef = er = -o.os, (3.23) 

as obtained earlier |1()2| ITKj . These values yield for the spin local charges 

= 0.647^ and = 0.747"', which explains the observed quenching of 
Ml strength. This quenching was one of the reasons for the long-standing 
problem of the " missing" Ml strength in heavy nuclei like ^°^Pb, which was 
greatly clarified by the polarized photon scattering experiments of Laszewski 
et al. |124j . (For details see the reviews ll2Uj and also ref. jSHI-) 

Such an approach, however, fails to explain the resonance widths that, 
especially in heavy nuclei, are caused by a coupling to more complex config- 
urations than those accounted for in the RPA. In other words, the RPA or 
TFFS calulations did not answer the important question of why there are so 
many 1+ levels with B(Ml)t< (2 - 3)/x^ observed in ^ospb [T1II[T21!. That 
was another part of the problem of the "missing" Ml strength. 

The existence of deviations from the RPA predictions was confirmed by 
the discovery of the low-lying Ml resonance in the magic nuclei ^^O (at 
about 15 MeV) 11221 and ^°Ca (at about 10 MeV) [Hi] because, accord- 
ing to the RPA, no Ml resonance with a similar energy can exist if the p 
shell in ^^O or the sd shell in ^°Ca is fully occupied. In order to illustrate 
this, we show in Fig. 3.10 the results of calculations of the Ml resonance 
in ^°Ca ini] performed within a simpler model of taking lplh®phonon 
configurations into account and without inclusion of the single-particle con- 
tinuum, as was briefly discussed in section 2.2. As was shown in refs. [221 ES], 
only the terms that correspond to the "backward-going" graphs — that is, the 
GSCphon graphs — are responsible for the effect. We found that inclusion of 
the GSCphon eliminated the prohibition of the existence of the low-lying Ml 
excitations in ^^O and ^'^Ca and, as one can see in Fig. 3.10, reasonably ex- 
plained at least the strong Ml excitations observed in ref. |125j . One should 
note, however, that the total Ml strength calculated is substantially larger 
than that observed |125j and that the calculations were performed using the 
Woods-Saxon form of the nuclear density in the ground state and with the 
phonons treated within the Bohr-Mottelson model. 

In Fig. 3.11 and Table 3.8 the results of the ETFFS calculations using 
Eqs. (3.22) and (3.23) for both isoscalar and isovector Ml excitations in 
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^°^Pb are presented. In spite of the above-mentioned quenching due to the 
spin local charges, it was necessary to include the complex lplh®phonon 
configurations as well as our GSCphon in order to obtain agreement with 
experiment. The width, which can be deduced from the curve, agrees with 
the experimental value of 1 MeV |124j (smearing parameter A = 100 keV). 
The Sji?j(Ml) 1 value of 11.57 /ig in the 6.3-8.7 Mev interval and the mean 
energy E = 7.7 MeV agree reasonably well with Eji?j(Ml) tex'p~ 15.6 /Xq and 
Eexp = 7.3 MeV |124j . A reasonable agreement for the isoscalar 1"*" level 
at E'exp = 5.85 MeV was obtained. A more detailed discussion of the Ml 
calculations, including those for ^^Ca and ^^Fe, can be found in refs. [201 IHH 

EH]. 

3.2.4 Effects of GSCphon 

We define the case of absence of GSCphon as that when the associated compo- 
nents of the full propagator A{uj), Eq. (2.54) or — in other words — when the 
quantities Q^-+\ Q^+'\ P^++\ Eq. (2.45) are equal to zero. In principle also 
the unlike components A*^ ^ give rise to GSCphon but it turns out that this 
effect is very small, as indicated in Tables 3.4, 3.5 and 3.6, where the third 
row includes the effects of the A^ -* only. These correlations are similar to 
the RPA GSC in the sense that they affect only the locations and values of 
the residues of the poles but do not change the EWSR. We will discuss in 
the following only GSCphon, which are created by the associated components 
of the full propagator and which gives some new effects, as were discussed in 
section 2.3.4. and demonstrated for the case of a "pure" GSC^'*"'^ effect in 
section 3.2.3. 

In all the calculations for ^''•^'^Ca and ^^Ni in the large energy intervals, 
which are given in Tables 3.4, 3.5 and 3.6, we obtained 97-102 percent of 
the corresponding EWSR for the RPA well as for that of the RPA + 

lplh®configurations (without the above-mentioned GSCphon)- Taking into 
account our GSCphon, however, increases the EWSR by 4-7 percent as a rule. 
This result is in accord with the result obtained analytically in ref. jnHl; al- 
though our model and the ESRPA used in that work differ greatly. The main 
differences are that we use complex configurations with collective phonons 
while in the ESRPA "pure" 2p2h configurations are used and, in contrast 
to the ESRPA, we do not use perturbation theory in the effective interac- 
tion. The complex configurations in the ground state give rise to an increase 
of < >p and, therefore, of the EWSR, and perhaps to changes in other 
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ground state characteristics. 

The role of GSCphon is especially noticeable in the low-lying energy region 
of the EO and E2 resonances in the nuclei under consideration. We have 
from Table 3.2 for the 5-12 MeV interval in ^°Ca 7 percent of our Sq and 
7.4 percent of Sj. The contribution of GSCphon to these figures is 3.6 percent 
and 4.8 percent, respectively — i.e, more than half. (For comparison, the 
contributions to this interval obtained within the RPA are 1.0 percent and 
0.4 percent). There is also additional low-lying strength due to GSCphon in 
^'^Ni (Fig. 3.4 ). There it is 4.8 percent of S[] in the 5.0-14.0 MeV interval 
as compared with 2.5 percent without GSCphon — i.e., about half. For higher 
energies the role of GSCphon is diminished; for example, for the 10-20.5 MeV 
interval in ^°Ca their contribution to the full E2 EWSR decreases it by about 
one-fourth. 

For Ml resonances, as can be seen from Table 3.8 and Fig. 3.12 for ^^Ni 
and ^^Ni and Fig. 3.11 for ^^'^Pb, the role of GSCphon is also noticeable, 
although probably not to the extent that it is for the electric resonances 
under discussion. This is in agreement with the results obtained for Ml 
resonances within the ESRPA [22]. Their role in the fine structure of the El 
strength in ^'^^Pb has been calculated [H^ and will be discussed briefly in the 
next subsection. 

3.2.5 The pygmy resonance and low-lying structures 

One of the impressive examples of the gross or fine structure in a broad re- 
gion near the nucleon separation energy is the so-called pygmy resonance. At 
present this subject is the subject of active discussion in connection with the 
general interest in nuclei with a large neutron excess because, phenomeno- 
logically, the pygmy resonance is described as the vibration of the neutron 
excess against the core with N = Z. Thus, properties of the pygmy reso- 
nance should strongly depend on the N/Z ratio |141j . 

Resonances of the same name were observed in (n, 7) spectra at neutron 
energies from 10 keV up to several MeV in many nuclei [53 I142j . They 
were described by means of an additional lorentzian in the radiative strength 
functions with the fitted parameters E ^ 2-6 MeV, P ~ 1.0-1.7 MeV and an 
integral strength of 0.1-1.0 percent of the classical sum rule jSSj- In the (7, n) 
cross sections in ^osp]^ there are well-known structures in the 7.6- 12.0 MeV 
interval [SH] and between 9.9 and 11.2 MeV |143j . which were partly manifest 
earlier in the photoabsorbtion cross section |105j too. 
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In different experiments for many nuclei with A = 58-208, resonance-like 
structures were observed as "the low-lying 52A~^/^ MeV dipole resonance" 
|54j . Its characteristics are F = 1.2-2.0 MeV and an integral strength of 
1.65-2.5 percent of the experimental integral strength in the region of the 
isovector El resonance. There are also other low-lying structures in the broad 
energy region near the nucleon separation energy; see, for example, ref. |H7] . 

In the broad energy region under consideration, one should observe at 
least E2, Ml and — mainly — El transitions. As yet, even for ^°^Pb, there are 
no consistent microscopic calculations with non-separable forces performed 
within the same calculational scheme for all three multipoles. It is clear that 
the usual RPA or QRPA calculations are not able to explain these structures. 
The calculations should take into account more complex configurations and 
the single-particle continuum for the energies above the threshold. Moreover, 
as calculations in ^osp]-, j-^^ye shown jSl], the role of GSCphon is also essential. 
The El photoabsorption cross section in the 4.5-7.5 MeV interval has been 
calculated there within the ETFFS, both with and without GSCphon (with 
smearing parameter A = 80 keV), and a very noticeable difference between 
these two cases was found. In particular, the known structure at 5.5 MeV was 
explained just by the GSCphon in these calculations. The experimental value 
of J]B{E1) = 1.338 e^-fm^ for six 1~ levels is in good agreement with the 
corresponding theoretical value of 1.304 e^-fm^. See also a discussion of this 
question in ref. [75] . 

The results of nuclear resonance fluorescence experiments should signif- 
icantly clarify the situation because they are able in principle to solve the 
long-standing problem of the identification of separate 1~,2"'" and l"*" states 
|144j . In connection with the numerous and successful results of these exper- 
iments, and with the interest in neutron-rich nuclei, the term "pygmy dipole 
resonance" (FDR) is used now for the El part of the "old" pygmy resonance. 

The new experiments in the neutron-rich oxygen isotopes ^^O, and 
^^O, which used electromagnetic excitation in heavy ion collisions at beam 
energies of about 600 MeV/nucleon, gave up to 10 percent of the classical El 
sum rule in the region up to 15 MeV |145j . which is in sharp contrast to the 
dipole response of stable nuclei where the experiments usually give about 1 
percent. It is clear that the role of the single-particle continuum is essential 
here but, as was discussed earlier, the role of complex configurations and, 
for non-magic nuclei, of a new (phonon) mechanism of pairing |146j should 
be important for such delicate properties of nuclei, especially of unstable 
ones. Because there are already many experimental data |1441 1148j for stable 



61 



nuclei, the comparison with the theory, which uses a non-separable nucleon- 
nucleon interaction and takes into account simultaneously the single-particle 
continuum and complex configurations, will enable the choice of a reasonable 
variant of the microscopic theory to apply to calculations in unstable nuclei. 

3.3 Giant resonances in unstable magic nuclei 

Because the ETFFS takes into account the single-particle continuum and 
uses the universal parameters of the effective interaction and local charges, 
which are the minimum necessary for calculating properties of unstable nuclei 
(see sections 1.3.2 and 1.3.3), this approach is apphcable to such calculations. 
On the other hand, the corresponding measurements could be a convincing 
verification of the universality of the ETFFS parameters. 

The results of the calculations of EO and E2 resonances in the unstable 
nucleus ^^Ni are presented in Figs. 3.4 and 3.5 and in Tables 3.5 and 3.6. 
The El calculations for this nucleus are given in Table 3.4 and in ref. jSH]- 
These results are largely similar to those for ^''Ca and ^®Ca discussed above. 

As discussed in section 1, the general interest in unstable nuclei is con- 
nected first of all with the neutron-rich nuclei. We have therefore calculated 
the isovector El [SI] and the isoscalar and isovector Ml resonance jSH ED] 
(see Fig. 3.12) in the very exotic nucleus ^^Ni. 

The calculations of the El resonance in refs. [SH ED] were performed us- 
ing our Woods-Saxon potential to calculate the density in the nuclear ground 
state in the interpolation formula (1.11), and thus with another value of the 
parameter /ex, but that is not critical — not, at least, for integral character- 
istics. It was found that the isovector El resonance in '^^Ni had a width of 
5.2 MeV. Its mean energy E = 15.2 MeV is much less than that obtained 
from the empirical formula E = 78A~^^'^ (= 18.3 MeV). In this case the in- 
fiuence of the lplh®phonon configurations is less than in other, "less exotic" 
nuclei: within the CTFFS we obtained E = 14.7 MeV and F = 4.8 MeV. 

As can be seen from Fig. 3.12, the Ml resonance in ''^Ni has an asymmetric 
shape with a width of about 1 MeV, whereas this resonance in ^^Ni and ^^^Sn 
has no fragmentation width. 

The Ml resonance in the neutron-deficient ^'^'^Sn nucleus was split into 
two major peaks, with Ei = 9.8 MeV and E2 = 10.5 MeV (note the smearing 
parameter A = 100 keV). This result was obtained without GSCphon, however 
in ref. |14()j they were taken into account in a simpler model and the result 
did not change. We see from Fig. 3.12 that these features for ^'^Ni and ^°°Sn 
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were caused by the inclusion of the complex configurations and for ^°°Sn also 
by the small value of the proton separation energy Bp = 2.91 MeV (for ^^Ni 
Bn = 5.98 MeV was used). 

4 Microscopic transition densities and the cal- 
culations of cross sections 

4.1 Comparison of microscopic and phenomenological 
transition densities 

The transition densities p^, which are necessary to describe the nuclear struc- 
ture in calculations of cross sections, are simply connected with our density 
matrix piif, E + tr]) determined in Eq. (2.57): 

pUr,^E) = Im / dEpL{r,E + zr]), (4.1) 

where T,B{EL) is the B{EL) value summed over the interval Ai?. 

The isovector El transition densities calculated in our approach for the 
large energy interval were obtained in ref. jHS]. There is no significant differ- 
ence between the continuum RPA and our full calculation. In all cases there 
are maxima of proton and neutron transition densities on the nuclear surface 
with approximately equal amplitudes and opposite signs. This corresponds 
to the isovector nature of these giant vibrations. The existence of the max- 
ima on the nuclear surface corresponds to the usual phenomenological models 
used, especially the Goldhaber- Teller model for ^°Ca and ^^Ca. We pointed 
out also that the transition densities defined in Eq. (4.1) are also large inside 
the nucleus, especially for the protons in ^^Ca and for the neutrons in ^^Ni, 
which may be important in the analysis of (e, e') experiments done at large 
momentum transfer. 

In the experiments, however, a smaller energy interval is usually investi- 
gated. In order to understand better the dependence of p{r) on the energy 
region, we have calculated it for the small energy intervals of 5 MeV [S3] and 
2 MeV (see Fig. 4.1). There is a considerable difference in the radial form of 
the transition densities as compared with the large interval, as well as with 
the different intervals under consideration. In this respect the microscopic 
ptr are very different from the phenomenological ones used in conventional 
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analyses, where they are taken to be independent of energy over the whole 
energy interval |157j . Such a strong dependence on the energy region needs to 
be checked experimentally. It is clear that this dependence may be important 
for the description of the gross structure of cross sections. 

In Fig. 4.1 we compare the phenomenological transition densities for the 
IS EO resonance in ^°Ca with the theoretical transition densities obtained 
for two 2 MeV intervals in the main part of the resonance. We see no- 
ticeable differences for the two intervals, both between the theoretical and 
phenomenological densities and between the theoretical densities themselves. 
The behavior of the microscopic and phenomenological densities near the nu- 
clear surface, to which the {a, a') cross sections are most sensitive, is very 
different in the two intervals considered. Inside the nucleus the behavior of 
the transition densities also differs strongly for the two intervals; for example, 
for the 17-19 MeV interval the neutron microscopic density has three nodes, 
whereas the IS EO phenomenological transition density has always just one 
node. 

Differential {a, a') cross sections calculated with the transition densities 
of Fig. 4.1 are shown in Fig. 4.2 (details of the calculations are given in sec- 
tion 4.2.2, below). The values of the ratio Mn/Mp of the nuclear transition 
momenta, which are taken to be 1 for ^°Ca in the phenomenological approach, 
are also given in Figs. 4.1 and 4.2. In the microscopic approach there is a no- 
ticeable difference in the Mn/Mp ratios, both with respect to themselves for 
the two intervals and between the phenomenological and microscopic densi- 
ties. The cross sections calculated with phenomenological densities taken at 
microscopic values of Mn/Mp are also shown. The various cross sections dif- 
fer most around zero degrees, where our theoretical cross sections are smaller 
by 5-8 percent compared to the phenomenological analysis. 

In order to demonstrate the role of the size of the energy interval, we 
compare two theoretical ^°Ca(Q!, a') differential cross sections for the IS EO 
giant resonance in Fig. 4.3. The dashed line is the result of using a transi- 
tion density averaged over the whole energy range, 12-30 MeV; the full line 
corresponds to the differential cross sections that were obtained by the pro- 
cedure discussed, performed for the 2 MeV bin and summed over the same 
energy range. As one can see, around zero degrees the two cross sections 
calculated with microscopic densities differ by nearly 25 percent. The dif- 
ference between the results with the phenomenological |157j and our 2 MeV 
microscopic transition densities is even larger. 

In addition to the differences in energy dependence, there are further im- 
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portant differences between the microscopic and phenomenological transition 
densities: (i) In the microscopic approach, the simple relations for the pro- 
ton and neutron components of the nuclear transition momenta (as discussed 
above) and transition densities, such as the ones with the ratio N jZ^ do not 
exist. This results in a specific energy dependence of the Coulomb-nuclear 
interference in total transition potentials, which may even change a destruc- 
tive interference into a constructive one. (ii) The multipole decompositions 
of the cross sections are also very different in the phenomenological approach 
as compared to the microscopic one. In the first case one tries to extract the 
corresponding multipole composition by fitting to various total and differen- 
tial experimental cross sections, whereas in the microscopic approach these 
multipole decompositions are determined within the theoretical model with 
known parameters. Indeed, our microscopic results deviate appreciably from 
the ones derived from phenomenological approaches. 

4.2 Microscopic analyses of inelastic electron and al- 
pha scattering experiments 

4.2.1 Electron scattering 

In order to compare with electron scattering experiments, in which neither 
of the isoscalar EO and E2 resonances are disentangled, it is necessary to 
sum the E2 and EO strengths in the same proportion that they enter into 
the longitudinal electron form factor [32]. Thus the quantity [dB{E2) / du + 
{2^ / 1^11)63 [E^]) / duj] should be calculated and compared with experiment. 

In our calculations of electron scattering in ^°Ca jHS] we obtained a good 
description of the experiments As can be seen in Fig. 4.4, the de- 

scription of the splitting of the isoscalar (E2+E0) strength into three peaks 
observed in ref. at about 12 MeV, 14 MeV and 17 MeV was obtained. 
The value of the E2+E0 theoretical isoscalar EWSR in the observed interval 
10.0-20.5 MeV, which is equal to [6581+(25/167r)3729] = 8436 e^-fm^MeV, 
agrees with the experimental value (7899 ± 1580) e^-fm^ MeV. 

Our calculations of the electromagnetic E2 strength in ^°Ca |HS1 gave a 
nearly uniform distribution in the 12.0-22.0 MeV interval. If we take the 
isoscalar part of it, however, the E2 strength shows more structure, which 
is washed out by the isovector contribution. In earlier experiments |108j the 
contribution of (23.5 ± 4.7)% of the isoscalar EO EWSR in the 10.5-15.7 MeV 
was determined by two independent methods: (a, ai) and (a, aiao) reactions. 
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A contribution of (30 ± 6)% in the 10.5-20.0 MeV interval just from (a, ai) 
reactions was also found. We have obtained good agreement with experiment 
for the first interval (20.7 percent of our Sq in the 10.5-15.5 MeV region) 
and somewhat more (50.2 percent of Sg) in the second interval. See also a 
discussion of the IS E2 and EO resonances in section 3.2.2. 

Thus the ETFFS describes reasonably well not only the summed strengths 
but also the gross structure of the E2 and EO strength extracted from electron 
and alpha scattering experiments with ^°Ca. A noticeable influence of our 
GSCp/ion was also found [S21- We shall discuss below some of the ETFFS 
results for cross sections. 



4.2.2 a particle scattering 

The calculations of the cross sections were carried out with the modified 
code DWUCK4. The modification was done in such a way that the mul- 
tipole transition potentials were constructed by single-folding the complex 
density-dependent Gaussian effective a-nucleon interaction |130j with our 
microscopic (Eq. (4.1)) following the prescription given in ref. |134j . The 
four parameters of the strength and range of the real and imaginary parts 
of the effective interaction were adjusted to the experimental cross sections 
of the low-lying 2^ and 3~ collective states [321. ■'■^ -^^S- 4.5 we show the 
results of these calculations for the 3^ and 2^ levels in ^^Ni. The agreement 
between the theory and experiment allows us to conclude that the choice of 
the above-mentioned four parameters is good, so that they have also been 
used in the cross section calculations in the 10-35 MeV interval. The pa- 
rameters of the optical model potential were taken from ref. |135j . As in the 
experimental analyses of refs. , contributions of the IS and IV El and 

IS EO, E2, E3 and E4 resonances, which were calculated within our approach 
in the region under consideration, have been taken into account. 

In our analysis of the (a, ai) cross sections the following procedure was 
used: The theoretical transition densities for each of the six EL resonances 
were analyzed over 5 MeV energy intervals for ^^Ni and these densities were 
used as input into the DWUCK4 code. For every energy interval the DWBA 
cross sections daL{EL,9)/dfl were calculated. The inelastic a spectrum was 
then obtained from the expression 

-XE,0)-)^ ,^.2 ,^ [El. 9). (4.2) 



dndE ' ' ' ^7rTL{E- ElY + (rL/2)2 dVt 



66 



Here the summation runs over the six multipolarities considered and the 
parameters El and Tl were obtained from a lorentzian fit to the calculated 
strength function Sl{E) defined in Eq. (3.1). 

For the cross section calculations in ^°Ca the same method was used, but 
it was necessary to use 2 MeV bins because the IS EO strength and the total 
spectra have much more structure. As in the case of ^^Ni, the contributions 
of the same six resonances were taken into account. 

In the calculations of the giant resonances we used a standard Woods- 
Saxon single-particle basis. The residual Landau-Migdal interaction given in 
Eqs. (1.10), (1.11) and (1.43), with the parameters of Eq. (1.42), was used. 
Additional information about numerical details can be found in section 3.1. 
In all these calculations the smearing parameter A = 500 keV was used. 

4.3 Calculation of the ^^Ni and ^'^Ca {a, ai) cross sec- 
tions; comparison with experiment. 

A general problem in nuclei with A < 90 is that the isoscalar monopole 
resonance is very broad and no longer concentrated in one single peak. In 
moving to lighter nuclei the role of the surface becomes more important than 
in heavy nuclei and the (a, «/) reaction is very sensitive to the nuclear surface. 
Moreover, for light and medium mass nuclei the single-particle continuum 
becomes very important. All of this has to be considered if theoretical models 
for giant resonances in lighter nuclei are to be developed. Therefore the 
conventional RPA approach to collective states has to be extended by the 
inclusion of surface modes and the single-particle continuum. 

Until recently there were several open questions concerning the exper- 
imental information about the IS EO resonance in nuclei with A < 90 
Pm I106j . In some cases the magnitude of the detected strength was much 
smaller than the EWSR limit. An important example in this connection is 
^^Ni, where originally only 32 percent of the EWSR was observed with inelas- 
tically scattered a-particles These authors used standard data analysis 
techniques with the same phenomenological transition densities for dif- 
ferent excitation energies. For comparison, the same type of experiment in 
^°Ca and a similar analysis by the same authors showed (92 it 15)% of the 
EWSR |128j . Such uncertainties might have serious consequences for nuclear 
matter compressibility and its applications to astrophysics. 

Our analysis within the ETFFS |13HI132'] of the same experimental data. 
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where microscopic transition densities were used instead, gave 71.4 percent of 
the IS EO EWSR in the observed 12.0-25.0 MeV interval, compared with the 
32 percent of ref. [22] • In these calculations the energy interval considered 
was divided into 5 MeV bins, for which the calculations were performed 
separately. The results of the calculations allowed the assumption |13H I132j 
that a part of the IS EO strength in ^^Ni might be hidden in the experimental 
background. 

In the meantime the 12-25 MeV energy range of the original experiment 
in ^^Ni has been extended [32]. Here the authors present results from the 
12.0-31.1 MeV excitation region. Compared with the original analysis [32], 
two new ingredients were included by the same authors in their new work 
|33j : (i) a nuclear reaction description following the methods of ref. |130j 
was used and (ii) the giant resonance peak obtained after subtraction of the 
continuum was divided into 15 energy intervals from 1 to 3 MeV, each of 
which was analyzed separately. 

In what follows we discuss the ETFFS analysis of the (a, a') experiments 
in ^^Ni and ^°Ca, including the newest ones [321 1129j for both nuclei, where 
the isoscalar giant resonance region was investigated. This analysis includes 
both of the necessary effects mentioned at the beginning of this subsection. 
Within this model we calculated the distribution of isoscalar strength of the 
isoscalar E0-E4 resonances and the corresponding transition densities. From 
these transition densities we obtained, in the standard way, (a, a') cross 
sections that we compare with the experimental data in ^°Ca in the observed 
8.0-29.0 MeV interval |128j and for the isotope ^^Ni in several energy regions: 
12-25 MeV [221, 12.0-31.1 MeV [32] and 12-35 MeV [1221 • 

4.3.1 ^^Ni results 

In Fig. 4.6 we compare our microscopic calculations with the newest experi- 
mental results, given in ref. [22|; ioi the IS EO and E2 strength distributions. 
For the observed 12.0-31.1 MeV interval we obtained a value of the mean 
energy (defined as mi/rrio) of the IS EO resonance of 19.9 MeV and find 
81.5 percent of the EWSR. The experimental data are (20.301o:m) MeV and 
(74l6°)%, respectively. For the "old" 12.0-25.0 MeV interval the new phe- 
nomenological analysis finds (58±6)% of the EWSR [22]; which is now much 
closer to the microscopic value of 71.4 percent |131j obtained for the same 
interval. Thus both analyses now give very similar results. 

As one can see from the lower part of Fig. 4.6, however, for the IS E2 
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strength the two approaches still arrive at quite different results. Our E2 
resonance mean energy value and the depletion of the IS E2 EWSR defined 
in the 10.0-20.0 MeV energy interval are 19.1 MeV and 47 percent, while the 
phenomenological analysis gives 16.1 MeV and (115±15)%. There is also a 
disagreement over the IS El and E3 strengths. Our results are presented in 
Fig. 4.7. (See also the discussion at the end of this subsection.) 

In Figs. 4.8 and 4.9 we compare the same experimental data with our 
theoretical results in a slightly different way. There we obtain in both cases 
good agreement between our theory and the data for the total cross sections. 

In Fig. 4.8 we show the role of the background and the contribution of 
the various giant resonances to the total cross section. We compare our the- 
oretical results (full line with dots) with the experimental cross section at 
6 = 1.08° (histogram) of ref. We obtained these data by subtracting 

an instrumental background (denoted by "Backgr. (2000)") from the orig- 
inal experimental spectrum. The theoretical cross section is the sum of six 
different multipoles, of which we only show EO and E2. The straight line in 
the lower part of the figure denoted by "Backgr. (1996)" corresponds to the 
analysis of ref. where the considered energy interval, 12-25 MeV, was 
smaller. In the original analysis with phenomenological transition densities, 
only 32 percent of the IS EO EWSR limit was found. With the improved — 
but still conventional — analysis, the authors of ref. |130j obtained about 50 
percent. 

In Fig. 4.8, we also compare the 12-25 Mev interval and the previous 
background subtraction with the new, extended interval and the experimen- 
tal data, where a different background has been subtracted. We find here 
that the total IS EO resonance cross section in the 12-25 MeV energy range is 
138.3 mb/sr and corresponds to 71.4 percent of EWSR. The area under the 
old background line in this region, which is included in our 71.4 percent of the 
EWSR, corresponds to 22 percent of the IS EO EWSR limit, or 42.6 mb/sr. 
This strength had been subtracted as a part of the background in ref. [22] 
and was therefore not included in the analysis of ref. |130j . If we extend the 
analysis to the larger 12-35 MeV interval our theoretical model predicts 89.6 
percent of the IS EO EWSR limit. For the comparison with experiment see 
ref. 

In Fig. 4.10 we compare our theoretical cross sections with the data at 
6 = 4.08°. We reconstructed the experimental cross section from Fig. 4 
and Fig. 1 of ref. [32] • Good agreement between the theory and experiment 
[221 is obtained. We also see that at this angle the IS E2 resonance and 
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higher multipoles contribute most to the cross section, whereas the monopole 
contribution is small. 

Thus one can conclude that the new experimental data for the IS EO 
resonance in ^^Ni are in good agreement with the microscopic calculations, 
which do not contain any fitting parameters for the nuclear structure part. 
The values of integral characteristics correspond now to the known experi- 
mental systematics. The disagreement between our microscopic results and 
the phenomenological analysis [32] for the IS E2 resonance only confirms our 
earlier conclusion |131j regarding the necessity of using microscopic transi- 
tion densities in the experimental analysis. The same is true for the IS El 
resonance, where the two approaches also lead to very different conclusions. 
The authors of ref. [22] obtained only 41 percent of the IS El EWSR and 
this strength was spread more or less uniformly from 12 MeV to 35 MeV. 
Our distribution of this strength is shown in Fig. 4.7. One can see that the 
distribution is not uniform; rather, it has a resonance structure. We obtained 
89 percent of the IS El EWSR and a value of 25.0 MeV for the mean energy 
in the interval under consideration. These figures are consistent with the 
results of ref. |136j for other nuclei. Our E3 IS strength is more uniformly 
distributed. 

4.3.2 ^°Ca results 

Our improved calculations of the IS EO resonance in ^'^Ca, presented in 
ref. 1^3, show that it has a more compact form than in our earlier calcu- 
lations |137j . but it remains strongly structured and spread out over a large 
energy interval: 65 percent of the EWSR is in the 11-23 MeV interval and 
106.7 percent is in the 5.0-45.0 MeV interval. Thus it is important that the 
large energy interval, 8.0-29.0 MeV, was studied in ref. |128j and that, in 
fact, most of the IS EO strength was found in their analysis. 

Our theoretical distributions of the IS EO and E2 strengths in ^°Ca are 
shown in Fig. 4.11. Compared with the corresponding results of the anal- 
ysis of ref. |128j . which were obtained from the difference between spec- 
tra taken at 9 = 1.1° and 6 = 2.4° (see Fig. 6 of ref. |128p . a reasonable 
agreement for the central part (10-23 MeV) of the IS EO resonance was 
obtained. The percentages of the IS EO EWSR in the four observed in- 
tervals 7.5-12.5 MeV, 12.5-22.5 MeV, 22.5-28.8 MeV and 7.5-28.8 MeV 
are 6.0(7.6±0.2)%, 60.0(50.0±1.4)%, 16.0(34.7±1.7)% and 81.6(92±2)%, re- 
spectively. (The quantities in parentheses are the results of ref. |128j with 
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statistical errors only.) The final result of the analysis in ref. |128j of the IS 
EO EWSR for between 8 and 29 MeV is (92±15)%, which agrees with 
our value of 81.6 percent. There is, however, noticeable disagreement in the 
low-lying and high-lying regions of the excitation spectrum. 

In ref. [HH] the authors find (33±4)% of the IS EO EWSR and (57±6)% of 
the IS E2 EWSR at a peak energy of 17.5±0.4 MeV. Our results for the IS EO 
and E2 resonance mean energies, obtained (as mi/mo) from averaging over 
the observed 8.0-29.0 MeV interval, are 17.2 MeV and 17.1 MeV, respectively. 
It is impossible to compare the experimental IS EO EWSR depletion value 
with our value of 81.6 percent, one of the reasons being that the experimental 
value of (33±4)% does not contain the continuum |128j . For the IS E2 
resonance, however, the difference between our value of 88.0 percent, obtained 
for the 8.0-29.0 MeV interval, and the 57 percent in ref. |128j is smaller than 
in the IS EO case, which may indicate that there is less E2 strength in the 
experimental continuum. 

A clearer comparison with experiment is shown in Fig. 4.12 for the double- 
differential cross section for ^°Ca. Our full calculations reproduce the exper- 
imental gross structure reasonably well but, again, we see differences at low 
and high energies. The general difference between the values of the theo- 
retical and experimental cross sections may be hinting that multipoles other 
than the IS EO contribute to that cross section. 

In order to understand better the role of complex configurations, we also 
show in Fig. 4.12 (a, a') cross sections obtained with the microscopic tran- 
sition density calculated within the continuum RPA only; that is, without 
inclusion of our complex IplhCSphonon configurations. It can easily be seen 
that there is a big difference between the two theoretical curves. We conclude 
that the gross structure of the IS EO resonance in ^°Ca is caused by the com- 
plex lplh®phonon configurations. A similar conclusion for the fine structure 
of the IS E2 resonance in ^''^Pb was reached in our earlier calculations [5^ 
and in section 3.2.2. 

As seen in section 4.3.1 and Figs. 4.8, 4.9 and 4.10, the calculations for 
^^Ni explained reasonably well not only the IS EO resonance in this nucleus 
but also the total spectra. We therefore calculated the total spectrum of 
the ^°Ca(a,a') reaction observed in ref. |128j . It has a detailed structure 
but, unfortunately, is given there only as counts. We obtained the following 
results, shown in Fig. 4.13: 

1. The ETFFS describes reasonably well the gross structure of total spec- 
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tra, except for the region below 9 MeV, which may be connected with 
the excitations of the a+^^Ar system. 

2. The difference between the two horizontal axes in Fig. 4.13 gives an 
instrumental background of the experiment under discussion. This dif- 
ference is about 7 mb/sr-MeV, so that the integrated cross section is 
about 7x25 = 175 mb/sr. (These numbers are, of course, very ap- 
proximate because they were obtained by imposing the experimental 
curve on the theoretical one). Integration of the theoretical curve gives 
422 mb/sr and of the experimental one 676 mb/sr. Thus we should 
compare (676-175) = 501 mb/sr (experiment) and 422 mb/sr (theory). 

3. As can be seen in Fig. 4.13, the experimental background shown by 
the dot-dashed line contains a noticeable contribution from giant reso- 
nances, including that from the IS EO. 

Thus, in the calculations presented in Fig. 4.12 and 4.13, a reasonably 
good description of the gross structure, both in the 0° cross section for the 
IS EO strength and in the total spectra in ^°Ca, has been obtained. One 
can, however, see some disagreement with the experiments around 19 MeV. 
Fortunately, the results of the new and improved experiments |129j have been 
published recently, and these should be analyzed within the ETFFS. 

Description of the improved EO cross section data 
in 40Ca(a, a') at = 240 MeV. 

According to the authors of ref . |129j , the new experimental results were 
obtained using an analysis technique that unambiguously identifies multi- 
pole strength, whereas the old spectrum subtraction technique |128j is very 
sensitive, particularly to experimental background. In ref. |128j a definite 
assignment could be made only for (33±4)% of the EO strength in the peak. 
What is also important is that the new method allows considerable extension 
of the observed energy interval to obtain thereby some additional EO strength 
in the high energy tail region, so that a new experimental background must be 
taken. In other words, the authors have extracted an additional EO strength 
from the old background, as was assumed in our work on ^^Ni beginning in 
1996 In the new ^^Ca measurements (97±11)% of the EWSR 

was observed in the 10-55 MeV energy interval, whereas the original interval 
was 4-27 MeV. This result agrees with our calculations: 81.6 percent in the 
(7.5-28.8) Mev interval jHZ] and 106.7 percent in the (5-45) MeV interval |HS]. 
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A possible reason for the above-mentioned disagreement with the ex- 
periment around 19 MeV could be that the number of low-lying phonons 
used — the 3~ and 5~) — is insufficient. We have therefore added a third one. 
Because this additional 2^-phonon at E = 3.90 MeV has another (probably 
2p2h) nature as compared with those used, we had to make further approx- 
imations for this specific phonon. We have therefore calculated, first of all, 
the distribution of the EWSR for the IS EO resonance with three phonons 
and compared it with experiment in Fig. 4.14. In Fig. 4.15 we compare the 
new |129j and old |128j experimental results for the EO cross section. One 
can see a noticeable difference between the experimental results of 1997 and 
2001. In Fig. 4.15 we show also our previous calculations within the ETFFS 
[H7j with two phonons, which were used earlier in the calculations for ^°Ca. 

With the inclusion of the third phonon, very good agreement with the 
new experimental results was obtained, both for the EWSR distribution 
(Fig. 4.14) and for the IS EO cross section (Fig. 4.16).^ 

We conclude that experiments of this type can shed some light on low- 
lying excitation spectra and, for a microscopic theory, can demonstrate the 
necessity of accounting for — and identification of — additional IplhCSphonon 
configurations, or even of inclusion of more complex configurations than the 
lplh®phonon. 

5 Conclusion 

In this review we presented a new microscopic many-body theory for the 
structure of closed shell nuclei. The extended theory of finite Fermi sys- 
tems (ETFFS) is based on the Landau-Migdal theory of finite Fermi systems 
(TFFS) and includes in a consistent way configurations beyond the Iplh 
level. A large part of this review is concerned with the application of this 
new approach to giant resonances in closed shell nuclei. As in the stan- 
dard TFFS j2] , we formulate the theory within the framework of many-body 
Green functions. As in the original approach, one makes an ansatz for the 
propagators in the equation for the response function that takes into account 
these higher configurations and one renormalizes the resulting equation in 
the standard way. After a long mathematical procedure one ends up with an 
equation that includes effective charges and an effective interaction which are 

■^Our estimates of the contribution of the fourth, i.e. , low-lying phonon at 3.35 
MeV, have shown that its contribution is small. 
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parameterized in the same way as in the standard theory. The corresponding 
parameters are universal for all the nuclei, investigated so far. The ETFFS 
extends the standard TFFS in the following directions: 

1. In addition to the Iplh configurations, it considers in a consistent way 
complex configurations of the lplh®phonon type. The corresponding 
formulation is a natural extension of the conventional RPA. 

2. It includes explicitly the single-particle continuum. For that reason one 
is able to obtain a microscopically determined envelope of the resonance 
without using a purely phenomenological smearing parameter. (For 
numerical convenience we still use a small smearing parameter of a few 
hundred keV.) This feature of the ETFSS is especially important for 
the nuclei with the nucleon separation energy near zero. 

3. The approach takes into account in a consistent way ground state cor- 
relations caused by the conventional RPA and by the more complex 
configurations. It is demonstrated that the latter ones (GSCp/jo^) are 
at least as important as the conventional RPA ground state correla- 
tions. 

Another important feature of the ETFFS is that the final equation explic- 
itly contains both the effective ph interaction and the quasiparticle-phonon 
interaction. This new interaction, however, does not introduce additional 
parameters. It is completely determined by the original ph interaction that 
enters into the conventional RPA equation. From this equation one obtains 
the structure of the phonons and the corresponding quasiparticle-phonon 
interaction. 

We have reviewed the main results obtained within the extended theory 
and discussed in detail the following physical results: 

1. Quantitative explanation of the widths of GMR. It is shown that the 
IplhCSphonon configurations give the most important contributions 
to the widths of the GMR. In addition, a large fraction of the ob- 
served gross and fine structure can be directly traced back to specific 
lplh®phonon configurations. Detailed results have been presented for 
^°Ca, "^^Ca, ^"^Ni, ^^Ni and ^"^Pb. These calculations give not only the 
solution to a long-standing problem, but they are also crucial for the 
analyses of the experimental data in medium mass nuclei. 
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2. Necessity of using microscopic transition densities: In order to extract 
the parameters of GMR from the experimental cross sections one needs 
the corresponding form factors. So far, in nearly all the analyses of 
hadron and electron scattering experiments, phenomenological transi- 
tion densities derived from some macroscopic models have been used. 
In heavy mass nuclei, where the GMR are relatively narrow and in gen- 
eral well separated from each other, such a procedure seems to work. 
For nuclei with A < 90, where the giant resonances show pronounced 
structure and are spread out over a large energy interval and where 
various multipoles overlap, such a procedure is no longer appropriate. 
First of all, the microscopic transition densities are energy dependent 
and vary over the ranges of the GMR in such nuclei. Secondly, as the 
resonances overlap, it seems difficult with the conventional method to 
extract in an unique way the width and the absolute strength of the 
individual GMR. Those effects were discussed in some detail for ^^Ni 
and ^"Ca. It was shown in particular that overlapping multipole reso- 
nances give rise to a very smooth angular distribution that looks like 
an experimental background. Such a misinterpretation gives rise to 
"missing strength," as occurred some years ago. The analysis of the 
corresponding experimental cross sections with microscopic transition 
densities derived from the ETFFS has clarified that point. 

3. It has been shown that the ground state correlations induced by com- 
plex lplh®phonon configurations give rise to a 5-7 percent contribu- 
tion to the EWSR. This additional strength violates the sum rule. One 
should note that this important question still lacks a completely satis- 
factory answer. 

In the next stage of the development of this new microscopic nuclear 
structure theory one has to include pairing. Such a generalization will be 
crucial for the application of the ETFFS to open shell nuclei. One may 
expect some new effects for the theory of nuclear pairing if one takes the 
quasiparticle-phonon interaction into account |146j . Such new effects could 
probably be observed in experiments with high experimental resolution. The 
first steps towards extending the ETFFS in this direction have already been 
taken |147j and the corresponding calculations are in progress. 

As mentioned earlier, the ETFFS and the original TFFS start with phe- 
nomenological single-particle potentials and parameterize the residual ph in- 
teraction; i.e., the theories are not self-consistent. For the applications of 
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these theories to nuclei far from stabihty the self-consistency, however, is 
crucial. Here one starts with, for example, a density functional or some 
relativistic Lagrangian, which provides a single-particle scheme and, simul- 
taneously, the residual interaction. Actually, the ETFFS is most appropriate 
for a self-consistent approach because it takes into account explicitly the ef- 
fects of the phonons on the single-particle energies and single-particle wave 
functions as well as in the residual interaction. In this respect the self- 
consistent ETFFS will go far beyond all existing theories for the structure 
of unstable nuclei. For closed shell nuclei, initial steps in this direction have 
already been taken |15Uj . We are convinced that for a quantitative under- 
standing of unstable nuclei one has to consider simultaneously non-separable 
universal forces, complex configurations, the single-particle continuum and — 
perhaps — the GSCphon, as discussed in section 1.3.2. We are also convinced 
that the new effects, such as the GSCphon (section 3.2.4) and the second 
(phonon) mechanism of pairing |146j . can manifest themselves in unstable 
nuclei in a clearer form than in the stable ones. The advent of modern 
radioactive beam accelerators and large gamma detector arrays will open 
possibilities to study these exciting new effects. 

The general outline for the construction of a theory with complex con- 
figurations used in the formulation of the nuclear ETFFS under discussion 
may be used, of course, for other Fermi systems — with and without super- 
conductivity, such as metallic clusters and quantum dots, where phenomena 
similar to zero sound have already been observed. 

Last but not least, all the formulas of the ETFFS can be straightforwardly 
generalized to the finite temperature case because the GF formalism used 
here has a natural generalization in the Matsubara temperature technique. 
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Table 1.1. Integral characteristics of the E2 IS and IV resonances 
in ^'^Ca,^°^Pb calculated within the Continuum TFFS 



Nucleus, 
resonance 


Interval, 
MeV 


EWSH 0/ 
EWSRth ' 


MeV 


i ,U 5 

MeV 


MeV 


MeV 


4"Ca, IS 


5-45 


100 


20.7 


21.6 


22.5 


~ 18 [Si! 


^^«Pb, IS 


0-45 


95 


6.1 


8.1 


10.2 


10.9 ± 0.3 jSH] 

10.4 inni 


208pb, IV 


0-45 


87 


16.6 


18.2 


19.3 


20.2 ± 0.5 jSZl 
22.6 ± 0.4 jSH] 
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Figure captions for section 1 

Fig. 1.1 Graphical representation of tlie Betlie-Salpeter equation for tlie 
response function in the ph channel. 

Fig. 1.2 The hadron strength function of the isoscalar E2 resonance in 
^°Ca calculated within the continuum TFFS (solid line). The dashed line 
gives the same quantity obtained without the ph interaction ("free response"). 

Fig. 1.3 The same as in Fig. 1.2 but for ^o^Pb. 

Fig. 1.4 The same as in Fig. 1.2 but for the isovector E2 resonance in 

208pjg 
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Figure captions for section 2 

Fig. 2.1 The energy-dependent mass operator and ph amplitude If^ 
in tlie g^ approximation. 

Fig. 2.2 Graphs corresponding to the propagator of the model IHH] 
including RPA (a) and lplh(8)phonon configurations (b,c,d). The wavy line 
and the circle denote the phonon Green function and the amplitude g, and 
the thin lines denote the Green function G. 

Fig. 2.3 Typical graphs of fourth order in the quasiparticle-phonon in- 
teraction amplitude g entering into the diagrammatic expansion in the time 
representation that correspond to the MCDD propagator used in the ETFFS. 
The notation is the same as in Fig. 2.2, but here the direction of the arrow 
on a fermion line denotes the particle or hole GF G. The dotted line denotes 
the time cut at fixed time. Only graphs d, e and f are explicitly included in 
the ETFFS. 

Fig. 2.4 Graphs of order g^ with self-energy insertions in the time repre- 
sentation that are included in the ETFFS and give a part of the GSGphon- 

Fig. 2.5 Typical "backward-going" graphs of order g^ due to the GSCphon- 
graphs a, b, c, and d are similar to the RPA GSC case, whereas graphs e, f, 
g and h are the GSCphon of new type that correspond to direct excitation of 
a complex configuration by an external field. 
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Figure captions for section 3 

Fig. 3.1 The El photoabsorption cross section for ^°Ca obtained witliin 
CTFFS (dasli-dotted), in tlie lpltL+lplli®plionon (daslied) and lpltL+lplli®phonon 
+ continuum (solid) approximation. The experimental data are taken from 
ref. [inSl. From [TBI. 

Fig. 3.2 Same as Fig. 3.1 but for '^^Ca. The experimental data are taken 
from ref. |lU4j . From [75] . 

Fig. 3.3 Same as Fig. 3.1, but for ^'^^Pb. The experimental data are taken 
from ref. |lU5j . From [75] . 

Fig. 3.4. The full EO and E2 electromagnetic strength functions calculated 
in the 5.0-45.0 MeV energy interval. The final results are shown by the solid 
curve. The RPA + lplh®phonon (without GSCphon) and the continuum 
RPA results are shown by the dashed and dotted curves, respectively. /^From 

Fig. 3.5. The full isoscalar and isovector EO and E2 strength functions cal- 
culated with taking into account the RPA and lplh®phonon configurations, 
the single-particle continuum and GSCphon- From [53] . 

Fig. 3.6 The EO isoscalar (full line) and isovector (dashed line) hadronic 
strength functions in 208p]^ calculated in the ETFFS (smearing parameter 
A = 250 keV). 

Fig. 3.7 The E2 hadronic strength function in ^osp]-, calculated within the 
ETFFS with smearing parameter A = 250 keV (upper panel) and 800 keV 
(lower panel) 

Fig. 3.8 The El photoabsorption cross section for ^°^Pb calculated in 
the RPA + lplh®phonon -|- continuum approximation with the smearing 
parameter A = 250 keV (heavy line) and A = 10 keV (thin line) compared 
with experiment. From [75] . 

Fig. 3.9. Upper: high resolution inelastic proton and electron scattering 
spectra in the region of the isoscalar E2 resonance in ^osp]^ [^ I86j . Lower: 
the isoscalar E2 strength function calculated in the RPA + lplh(8)phonon + 
continuum approximation. See text for details. 

Fig. 3.10 Ml in ^°Ca calculated in the RPA + lplh(8)phonon approxi- 
mation [HSlinHj- The dashed lines represent the experiment (s) [125j . From 

m- 

Fig. 3.11. Upper: the strength function for the Ml excitation in 208p|-j 
in the CRPA (dotted line) and in the RPA -|- lplh(8)phonon + continuum 
without (dashed line) and with GSCp^o^ [HUj (full line). Lower: experimental 
data [1211 . 

Fig. 3.12. The Ml strength function for unstable nuclei calculated within 
the CTFFS (dotted curve) and the ETFFS without (gs-) and with (gs+) 
GSPp;,^,. From [50] 
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Figure captions for section 4. 

Fig. 4.1. The microscopic (dotted line) and plienomenological |157j (solid 
line) IS EO transition densities in ^°Ca, calculated for two different energy in- 
tervals of 2 MeV each. Results for proton and neutron microscopic transition 
densities are also given (denoted by triangles). The values of the correspond- 
ing ratio Mn/Mp of the nuclear transition moments are given in the legend. 
From [HZI. 

Fig. 4.2 Angular distributions for the IS EO resonance in the same 
intervals as in Fig. 1, calculated with phenomenological |157j and microscopic 
transition densities. For comparison, the results with the phenomenological 
transition density, obtained with the microscopic value of the ratio Mn/Mp 
of the nuclear transition moments (dashed line) are also shown. From [H7j . 

Fig. 4.3. Angular distributions for the IS EO resonance in ^''Ca in the 12- 
30 MeV interval, calculated using two different methods: the division of the 
interval into 9 bins of 2 MeV each (solid line) and the interpretation of the 
12-30 MeV interval as one bin (dashed line), respectively. For comparison, 
the results obtained with the phenomenological transition density |157j are 
also given (solid curve with dots). From [H7] . 

Fig. 4.4 (a) The experimental data |lllj for the strength function vs. the 
theoretical one for ^°Ca (see text), (b) The theoretical electromagnetic E2 
(dashed line) and EO (dotted line) strength functions. From |83j . 

Fig. 4.5 The experimental [221 and theoretical a particle scattering dif- 
ferential cross sections for the 3J~ (below) and 2^ (above) levels in ^^Ni. 

Fig. 4.6. Distribution of the IS EO and E2 strengths in ^^Ni. The experi- 
mental data are taken from [33] • From |87j . 

Fig. 4.7. Distribution of the IS El and E3 strengths in ^^Ni as calculated 
within our theoretical model. From |87j . 

Fig. 4.8. Cross sections of ^^Ni(a,a') at E„ = 240 MeV and 9 = 1.08°. 
The experimental data (histogram) were taken from Ref. [221; where an in- 
strumental background has been subtracted. The solid curve with dots shows 
the calculated total (summed) cross sections, the dashed line ("Rest(EL)") 
corresponds to the sum of the IS and IV El, and IS E3 and E4 multipoles. 
The shaded area shows an additional IS EO strength which has been sub- 
tracted in the previous experiments as background [22] • This area corre- 
sponds to 22 percent of the IS EO EWSR, see text. From jHTj. 

Fig. 4.9. Cross sections of ^^m{a,a') at = 240 MeV and 9 = 1.08°. 
The experimental data (histogram), including the instrumental background, 
are taken from ref. jHSl- The solid curve with dots gives the calculated total 
(summed) cross sections. In the lower part of the picture the components 
of the total cross sections are shown without the background. In particular, 
the dashed line gives the IS EO contribution. From [S7j . 

Fig. 4.10. Same as in Fig. 4.9 but for 9 = 4.08°. The experimental data 
were obtained by using the results of ref. [S2] (see text). From |87j . 

Fig. 4.11 Distribution of the IS EO and E2 strengths in ^°Ca (theory). 
From [HZI. 



92 



Fig. 4.12. The 0° cross section for the IS EO strength in ^°Ca(a,a') at 
Ea = 240 MeV, calculated with (solid line) and without (dotted line) taking 
into account complex lplh®phonon configurations. The experimental data 
(histogram) are taken from ref. |128j . One can see that the gross structure 
of the IS EO resonance is caused by complex lplh®phonon configurations. 
From [H7I. 

Fig. 4.13 Cross sections of ^°Ca(a,a') at E„ = 240 MeV and 6 = 1.1°. 
The experimental data dotted line) and the background (dot-dashed line) are 
taken from ref. |128j . Because in Fig. 2 of ref. |128j there are only counts, 
the experimental histogram was imposed on the theoretical curve in such a 
way that the maxima of the two curves coincided. In this way it was possible 
to estimate roughly the contribution of the instrumental background as the 
difference between two horizontal axes (see text). 

Fig. 4.14 Upper: distribution of the IS EO EWSR calculated with 2 and 
3 low-lying phonons. Lower: comparison of the 3-phonon calculation with 
experiment |129j . 

Fig. 4.15 EO cross sections of ^°Ca(a, a') at E„ = 240 MeV and 6 = 1.1°. 
The experimental data 2001, 1997 and the theoretical results with 2 phonons 
are taken from refs. |1291 1128] and [HZ] , respectively. 

Fig. 4.16 Comparison of the improved experimental data |129j for cross 
sections and the theoretical calculations with 3 phonons. 
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Appendix. 

1. Generalized propagator in the method of chrono- 
logical decoupling of diagrams 

The equation for the density matrix /?i2(a;), Eq. (2.19), in the A representation 
is given by (the tilde on ex and fix has been omitted in the equations that 
follow) 

5pi2{uj) = 5p?2(^) - Yl ^12,34^34,56^P56(t^) 
56,34 

^P'M = - E^12,34(c^)(e,V^°)43 (A.l) 
34 

The generalized propagator A has the form 

^12,34 (^) = [^15^26 + <5^2756 (^)] ^56778 (^) [^73^4 + Q78!34(^)] + Pu%i^) ■ 

56,78 

(A.2) 

For the ph-ph part A^^ ^8 of the propagator the following equation must 
be solved 

^r2734(^^) = ^12,34(t^) - J2 ^12,56(t^)*56,78(t^)^78734(^^) (^-3) 

56,78 

where the propagator 

^12,34(t^j = ^13^24 

u - ei2 

ei2 = ei - €2, Hi = l{ei < eF),ni = 0(ei > ep), (A.4) 
is the RPA propagator and 

$12,34 = $12,34 + ^T27i (A.5) 
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-1) 



J+Ja 



5 a; - esr - (1 - 2ni)ug 



{ J, 


jl 






J4 


J3 J 



(l-ril)(l-n3) 

UJ - €32- U; + €41 + 

(1 - ni){l - n-i) _ n^n?. 

U — ei4 — Uq UJ + 623 + 



+ 



■1) 



J+Ja 



X ^ (^6 - n2)[uJ - 612 - £64 - (1 - '2ni)Uq] 

\ [662 + (1 - 2n,)ug\[e,, + (1 - 2n,)ug\ ^^4,6+ 
c ^ (ni - r25)[a; - 612 - 635 - (1 - 2^1)0;^] 

^^24 V [615 + (1 - 2ni)^,] [635 + (1 - 2ni)^,] ^13,5 + 

(1 - 711)713 (1 - n3)ni \ 



[ J, 


jl 


M 




J4 


J3 J 



+ 



(1 - 721)^3 (1 - n3)rai \1 -1 

624 + ^5 J' 



£31 + y 

(A.6) 



^S™4 = K -ni)[(l -ni)(l -^3) +721^3] [(1 -n2)(l -^4) +n2n4]x 

S /^li, 5/^24,6 ^ 

91, 92, 56 

(1 -72i)(l -?^5)(1 -^e) +nin5nQ 

662 + (1 - 2ni)t^,J[664 + (1 - 27li)uJg^\[uJ - 652 - 654 " (1 " 27lx){ijJq^ + C^gJJ 

(1 - r22)(l - n^){X - ne) + 722^5^6 

[615 + (1 - 2ni)t^gJ[635 + (1 - 27ll)uJg^][uJ - 615 - 636 - (1 - 27ll){ujg, + iOg,,)] 

[(1 - 7l2){l - 7l^)7lQ + 722^5(1 - rig)] 



[615 + (1 - 2ni)t^qJ[635 + (1 - 2rai)t^gJ[662 + (1 - 27li)ujg^] 
- ei2 - £34 + £56 - (1 - 27li){ujg^ +^92)] 



X 



[664 + (1 - 2ni)u;, 



92 J 



■}• 



(A.7) 
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The quantities and are given by 

Qi2^34 = [(1 -^^2) +nin2][{l -n3)ni + n-i{l -n4)]x 



q I 6 

+ 



1 — ni 



+ 



ni 



e24(e64 + t^g) ' e62+t^g^e24 



ni 



e42(e46+l^g) <^26 + '^g ^ <^42 iO - eiQ - UJg 



^^6/^24,6 



-^24 E 

5 



ni 



1 - 1 



71^(3 



?(-)* 

13,5 " 



1 — rii 



Til 



£31 (^53 + t^g) ^51 +^^^^^31 — 652 — LOq 



(l-^5)/5lT5 



+ 



-l)^'^"^«2(l-^l) 



(631 + CJg) (u; - 632 - UJq) (e42 + UJq) {uj + €41 + Uq) 
(1 1^3 



+ 



{e2i + ujq){uj -eii, - UJq) (ei3 + a;g)(a; + e23 + u;q) 

(A.8) 



Ql2'l34 — Q34.I2 



-fl2,34(^) 



E{5i3E 

n ^ a 



[(1 - r2i)(l - 722) + nira2] [(1 - n3)(l - ^4) + 7237^4] x 

(rii - n6)/3|4,6 

't [^26 + (1 - '^ni)ujq] [646 + (1 " 2ni)u;g] [t^ - ei6 + (1 - 2ni)t^5] 



5 [ei5 + (1 - 2ni)cu,] [e35 + (1 - 2ni)cc.g] [uo + €25 + (1 - 2ni)cc.,] 

\ Jq 3A J3 



+ 



ni{l - n^) 



1 



(e3i + UJq){e42 + CUq) + ^41 + UJ - 632 - UJq 



+ 



^3(1 - rii) 



1 



(ei3 + UJq) (624 + UJq) \UJ + €23 + UJq UJ - 614 



(A.9) 
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Here the following notations have been used: 

(A. 10) 



Pl2,3 — 2ji + 1 ^13y23) Pl2,3 — 2ji + 1 ^31532) 



A\3 = (l-^3)/5S+n3/3?y. 

Into all these formulae the reduced matrix elements =< 1||5'^||2 > 
(g = JLS) of the phonon creation amplitude at J = L. 

gir,cT)= Y: 9jLsir)Tj[sin,cT) (A.ll) 

JLSM 

enter, where n = r/r and 

Tjlsin, ct) = J2iLf^Su\JM)YL,in)as. (A.12) 

IJ,,U 

with u = ±1,0, aot. = 1 and = a^, (Tj, being cyclic components of spin 
matrices. In the spherical single-particle basis they are determined by 

< m,hj,n,\gjLsT^^s\^2l2j2n2 >= (-lyn+J+rm ^ 

X ( ^ <1||/^^||2>, (A.13) 

y -mi M 7712 J M ) V / 

where 

< 1||5^''^||2 >=< liJi\\TjLs\\l2j2 > J Rhnni9JLs{r)Ri,j,nydr. (A.14) 

The reduced matrix element of the spherical tensor operator Tjls has the 
form 

< hjl\\TjLs\\l2j2 > = 



= 1[1 + (- l^L+h+hT^^ ^/ (2J+l)(2L+l)(2ji+l)(2j;+l)" ,, 

^ ^ ^ \ \ (t. \ ( -[]J+L+s^\ I sjs+i) ( J L S \l( ]i j2 J 

j^^^^^ ' ""^'y [ 1 -1 7^1/2 -1/2 

(A.15) 

where Xk = {k - jk){'2jk + 1). 



2. The refinement procedure 

Eq. (3.18) for the refined single-particle energies eA has the form 

^~i = ^i-E , , (A.16) 

ei - 63 - (1 - 2773)UJg 

where index 1 stands for the set of the single-particle quantum numbers and 
ei are the phenonomenological single-particle energies. 
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